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PREFACE 



Industrial Biotechnology is an interdisciplinary topic to which 
tools of modern biotechnology are applied for finding proper 
proportion of raw mix of chemicals, determination of set points, 
finding the flow rates etc., This study is significant as it results 
in better economy, quality product and control of pollution. The 
authors in this book have given only methods of industrial 
biotechnology mainly to help researchers, students and chemical 
engineers. Since biotechnology concerns practical and diverse 
applications including production of new drugs, clearing up 
pollution etc. we have in this book given methods to control 
pollution in chemical industries as it has become a great health 
threat in India. In some cases, the damage due to environmental 
pollution outweighs the benefits of the product. 

This book has six chapters. First chapter gives a brief 
description of biotechnology. Second chapter deals will proper 
proportion of mix of raw materials in cement industries to 
minimize pollution using fuzzy control theory. Chapter three 
gives the method of determination of temperature set point for 
crude oil in oil refineries. Chapter four studies the flow rates in 
chemical industries using fuzzy neutral networks. Chapter five 
gives the method of minimization of waste gas flow in chemical 
industries using fuzzy linear programming. The final chapter 
suggests when in these studies indeterminancy is an attribute or 
concept involved, the notion of neutrosophic methods can be 
adopted. The authors feel that the reader should be well versed 
with fuzzy models like neural networks, fuzzy relational 
equations, fuzzy control theory, fuzzy linear programming and 
neutrosophic fuzzy models like NRE together with a knowledge 
of the technical functioning of chemical industries. 

The authors are deeply indebted to Dr. Kandasamy, Kama 
and Meena for their sustained cooperation. 

W.B.VASANTHA KANDASAMY 
FLORENTIN SMARANDACHE 



Chapter One 



Introduction 



In keeping with the definition that "biotechnology is really no 
more than a name given to a set of techniques and processes", 
the authors apply some set of fuzzy techniques to chemical 
industry problems such as finding the proper proportion of raw 
mix to control pollution, to study flow rates, to find out the 
better quality of products. We use fuzzy control theory, fuzzy 
neural networks, fuzzy relational equations, genetic algorithms 
to these problems for solutions. 

When the solution to the problem can have certain concepts 
or attributes as indeterminate, the only model that can tackle 
such a situation is the neutrosophic model. The authors have 
also used these models in this book to study the use of 
biotechnology in chemical industries. 

The new biotechnology revolution began in the 1970s and 
early 1980s when scientists learned to precisely alter the genetic 
constitution of living organisms by processes out with 
traditional breeding practices. This "genetic engineering" has 
had a profound impact on almost all areas of traditional 
biotechnology and further permitted breakthroughs in medicine 
and agriculture, in particular those that would be impossible by 
traditional breeding approaches. 



There are evidences to show that historically biotechnology 
was an art rather than a science, exemplified in the manufacture 
of wines, beers, cheeses etc. It is well comprehended by one and 
all that biotechnology is highly multi disciplinary, it has its 
foundations in many fields including biology, microbiology, 
biochemistry, molecular biology, genetics, chemistry and 
chemical and process engineering. It is further asserted that 
biotechnology will be the major technology of the twenty first 
century. 

The newly acquired biological knowledge has already made 
very important contributions to health and welfare of human 
kind. 

Biotechnology is not by itself a product or range of 
products; it should be regarded as a range of enabling 
technologies that will find significant application in many 
industrial sectors. 

Traditional biotechnology has established a huge and 
expanding world market and in monetary terms, represents a 
major part of all biotechnology financial profits. 'New' aspects 
of biotechnology founded in recent advances in molecular 
biology genetic engineering and fermentation process 
technology are now increasingly finding wide industrial 
application. 

In many ways, biotechnology is a series of embryonic 
technologies and will require much skilful control of its 
development but the potentials are vast and diverse and 
undoubtedly will play an increasingly important part in many 
future industrial processes. 

It is no doubt an interaction between biology and 
engineering. The developments of biotechnology are proceeding 
at a speed similar to that of micro-electronics in the mid 1970s. 
Although the analogy is tempting any expectations that 
biotechnology will develop commercially at the same 
spectacular rate should be tempered with considerable caution. 
While the potential of new biotechnology cannot be doubted a 
meaningful commercial realization is now slowly occurring and 
will accelerate as we approach the end of the century. New 
biotechnology will have a considerable impact across all 
industrial uses of the life sciences. In each case the relative 



merits of competing means of production will influence the 
economics of a biotechnological route. There is no doubt that 
biotechnology will undoubtedly have great benefits in the long 
term in all sectors. The growth in awareness of modern 
biotechnology parallels the serious worldwide changes in the 
economic climate arising from the escalation of oil prices since 
1973. 

Biotechnology has been considered as one important means 
of restimulating the economy whether on a local, regional 
national or even global basis using new biotechnological 
methods and new raw materials. Much of modern biotechnology 
has been developed and utilized by large companies and 
corporations. 

However many small and medium sized companies are 
realizing that biotechnology is not a science of the future but 
provides real benefits to their industry today. In many industries 
traditional technology can produce compounds causing 
environmental damage whereas biotechnology methods can 
offer a green alternative promoting a positive public image and 
also avoiding new environmental penalties. 

Biotechnology is high technology par excellence. Science 
has defined the world in which we live and biotechnology in 
particular will become an essential and accepted activity of our 
culture. Biotechnology offers a great deal of hope for solving 
many of the problems our world faces!. As stated in the 
Advisory Committee on Science and Technology Report 
Developments in Biotechnology, public perception of 
biotechnology will have a major influence on the rate and 
direction of developments and there is growing concern about 
genetically modified products. Associated with genetic 
manipulation are diverse question of safety, ethics and 
welfare. 

Public debate is essential for new biotechnology to grow up 
and undoubtedly for the foreseeable future, biotechnology will 
be under scrutiny. We have only given a description of the 
biotechnology and the new biotechnology. We have highly 
restricted ourselves from the technical or scientific analysis of 
the biotechnologies as even in the countries like USA only less 
than 10% of the population are scientifically literate, so the 



authors have only described it non-abstractly and in fact we are 
not in anyway concerned to debate or comment upon it as we 
acknowledge the deep and dramatic change the world is facing 
due to biotechnology and new biotechnology. 

For more of these particulars please refer [1, 2, 13, 15, 17]. 
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Chapter Two 



Biotechnology in Chemical 
Industries 



The chemical industries have become a great threat in India. For 
the problems they cause on environmental pollution is much 
more than the benefit derived by their product. Some of them 
damage other living organisms like fishes, plants and animals; 
some cause health hazards to people living around the industries 
like respiratory ailments, skin problems and damage to nervous 
systems. So we have chosen to illustrate the minimization of 
pollution by CKD in cement Industries. Most of these problems 
can be controlled provided one takes the proper proportion of 
the mix of raw materials, which would minimize the pollution. 

Cement kiln dust (CKD) emits nitrogen, carbon etc., that are 
pollutants of the atmosphere and the waste dust affects the 
smooth kiln operation of the cement industry system and it 
reduces the production of clinker quality. Hence the 
minimization of waste CKD in kiln is an important one in the 
cement industry. The control of the waste CKD in a kiln is an 
uncertainty. Researchers approach this problem by 
mathematical methods and try to account the waste CKD in a 
cement kiln. But, most of their methods do not properly yield 
results about the waste CKD in kiln. Further, the control of the 
waste CKD in kiln is a major problem for this alone can lead to 
the minimization of atmospheric pollution by the cement 
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industry. So in this chapter we minimize the waste CKD in kiln 
and account for the waste CKD in kiln using fuzzy control 
theory and fuzzy neural networks. 

In this chapter fuzzy control theory (FCT) is used to study 
the cement kiln dust (CKD) problem in cement industries. 
Using fuzzy control method this chapter tries to minimize the 
cement kiln dust in cement industries. Cement industries of our 
country happens to be one of the major contributors of dust. The 
dust arising in various processing units of a cement plant varies 
in composition. In 1990 the national average was 9 tons of CKD 
generated for every 100 tons of clinker production. The control 
of cement kiln dust is a very important issue, because of the 
following reasons : 1. CKD emits nitrogen, carbon etc., which 
are pollutants of the atmosphere, 2. The waste dust affects the 
smooth kiln operation of the cement industry system and it 
reduces the production of clinker quality. The following creates 
mainly this waste dust in three ways in cement industries : (a) 
Cement kiln dust when not collected in time and returned into 
the kiln, cause air pollution, (b) Process instability and 
unscheduled kiln shutdowns and (c) Mixing of raw materials. 

The data obtained from Graft R. Kessler [12] is used in this 
chapter to test the result. After using the data from Kessler [12] 
this chapter tries to minimize the CKD in cement factory. The 
minimization of CKD plays a vital role in the control of 
pollution in the atmosphere. 

W.Kreft [21] used the interruption of material cycles 
method for taking account and further utilization of the waste 
dust in the cement factory. But this method does not properly 
account the waste CKD. Kesslar [12] has used volatile analysis 
to reduce CKD. In the volatile analysis method the alkali ratio is 
used to indicate the waste amount of CKD in clinker. 

Kesslar [12] classifies the raw data under investigation in 
four ways : 

I. Monitor and control of the system 
II. Burning zone and fuel combustion improvements 

III. CKD reprocessing 

IV. Find the mix of raw materials in proper proportion. 

The ratio of alkali should be lying between 0.5 to 1.5 in 
Kiln load material. But in this method the CKD was 



12 



approximately estimated up to 40%. He has not exactly 
mentioned the percentage of CKD according to the alkali ratio 
in an online process. So this method has affected largely the kiln 
system. 

In this chapter, in order to account for the waste CKD, the 
variables are expressed in terms of membership grades. This 
chapter considers all the four ways of waste CKD mentioned by 
Kesslar [12] and converts it into a fuzzy control model. This 
chapter consists of five sections. In section 1 we describe the 
cement kiln system and the nature of chemical waste dust which 
pollutes the atmosphere. In section 2 we adopt the fuzzy control 
theory to monitor and control the system and give suggestion 
for the improvement of burning and combustion zone. Section 3 
deals with the determination of gas volume set point and 
temperature set point for CKD reprocessing which is vital for 
the determination of percentage of net CKD. The amount of 
waste dust depends largely on the mix of raw materials in 
proper proportion of raw material mix is shown in section 4. 
The final section deals with results and conclusion obtained 
from our study. 



2.1 Description of waste CKD in cement kiln 

The data available from any cement industry is used as the 
information and also as the knowledge about the problem. This 
serves as the past experience for our study for adapting the 
fuzzy control theory in this section. This chapter analysis the 
data via membership functions of fuzzy control method and 
minimizes the waste CKD in cement industries. Since the 
cement industry, emits the cement kiln dusts into the 
atmosphere, this waste dust pollutes the atmosphere. 

This analysis not only estimates the cement kiln dust in 
cement industries but also gives condition to minimize the 
waste CKD so that the industry will get maximum profit by 
minimizing the waste CKD in cement industry. 

CKD is particulate matter that is collected from kiln exhaust 
gases and consist of entrained particles of clinker, raw materials 
and partially calcined raw materials. The present pollution in 
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environment is generated by CKD along with potential future 
liabilities of stored dust and this should make CKD reduction a 
high priority. Here we calculate and minimize the net CKD in 
kiln system. This chapter tackles the problem of minimizing 
waste CKD in kiln system in four stages. At the first stage we 
monitor and control the system. In the second stage we adopt 
time-to- time improved techniques in burning zone and 
combustion. At the third stage CKD reprocessing is carried out 
and in the fourth stage we optimize the mix of raw materials in 
proper proportion using fuzzy neutral network. The above stage- 
by-stage process is shown in the following figure 2.1.1. Fuzzy 
control theory and fuzzy neutral network (FNN) is used in this 
chapter for the above - described method to minimize the CKD 
in kiln system. 



Step 4: Optimize 
and mix the raw 
material in proper 
proportion 



Step 3: CKD 
Reprocessing 




Step 1: Monitor 
and control of the 
system 



Step 2: Burning 
zone and 
combustion 
improvement 



FIGURE 1: CKD Reduction using fuzzy control 



The fuzzy controller is composed of linguistic control rule, 
which are conditional linguistic statements of the relationship 
between inputs and outputs. One of the attractive properties of 
fuzzy controller is its ability to emulate the behaviour of a 
human operator. Another important characteristic of a fuzzy 
controller is its applicability to systems with model uncertainty 
or even to unknown model systems. The use of fuzzy control 
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applications has expanded at an increasing rate in recent years. 
In this chapter we use fuzzy control to monitor waste dust in 
cement kiln system and CKD reprocessing. The fuzzy control in 
kiln system is described in the figure 2.1.2. We use fuzzy neural 
network method and tries to find a proper proportion of material 
mix in cement industries. 

The authors aim to achieve a desired level of lime saturation 
factor (LSF), silica modulus (SM) and alumina modulus (AM) 
of the raw mix, to produce a particular quality of the cement by 
controlling the mix proportions of the raw materials. To achieve 
an appropriate raw mix proportion is very difficult, due to the 
inconsistency in the chemical composition ratio given for the 
raw materials. 

Fuzzy neural network model is used to obtain a desired 
quality of clinker. The raw mix as per the norms of cement 
industries should maintain the ranges like LSF 1.02 to 1.08, SM 
2.35 to 2.55 and AM 0.95 to 1.25, which are the key factors for 
the burnability of clinker to obtain a good quality of cement. 
Fuzzy control theory method is used to minimize waste cement 
kiln dust. Fuzzy control theory allows varying degrees of set 
membership based on a membership function defined over a 
range of values. The membership function usually varies from 
to 1. 
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FIGURE 2: Fuzzy control in kiln system 
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2.2 Monitoring and control of the system using FCT 
and improvement of burning zone and combustion 

Monitoring and control of the system is the most effective 
method towards CKD reduction in environment. CKD consists 
mainly of raw materials, which contain volatile compounds, 
therefore, tracking and control of the volatile compounds 
throughout the system often allows for the minimal CKD. The 
initial step in our plan towards CKD reduction is to identify the 
amount of the CKD. Here the indirect weighing method is 
applied to identify the amount of the CKD. Calculating 
sulphur/alkali ratio is a good indication of a possible imbalance. 
This ratio is calculated as the molar ratio of S0 3 /(K 2 0)+Na 2 0) 
in kiln load material. 



CKD VOLATILE ANALYSIS 



Volatile 


Molecular Weight 


Na 2 


62 


K 2 


94.2 


S0 3 


80 



Ratio of alkali = S0 3 /K 2 + Na 2 = 80/156.2 = 0.512 

This ratio should be between the values 0.5 to 1.5 in Kiln 
load material. The industry knows upto 40% of CKD exits, 
when the alkali ratio is between the values 0.5 to 1.5. But they 
cannot say exactly how much percentage of CKD waste comes 
from kiln by using the ratio of alkali in the online process. If 
industry knows this correct percentage of CKD in the online 
process, they can change some condition in the kiln and thus 
reduce the CKD in the online process. We adopt fuzzy control 
to estimate the percentage of CKD by using the ratio of alkali. 
The alkali ratio, kiln load material in tons and percentage of 
CKD are measured from the past happening process in kiln on a 
scale from 0.5 to 1.5, 5 to 25 tons and to 40% respectively. 
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That is we assign the sulphur/alkali ratio shortly termed as alkali 
ratio, alkali ratio to be approximately low (L) when its value is 
0.5, medium (M) when its value is 1 high (H) when its value is 
1.5. In a similar way we give kiln load material = {5 tons [first 
stage (FS)], 15 tons [second stage (SS)] and 25 tons [third stage 
(TS)]}. Percentage of CKD = {0 [very less (VL)], 10 [less (L)], 
20 [medium (M)], 30 [high (H)] and 40 [very high (VH)]}. ('=' 
Denotes approximately equal). The terms of these parameters 
are presented in figures 2.2.1, 2.2.2 and 2.2.3. 




0.5 



Legend 
MSG - Membership grade 
L-low, M-medium, H-high 



1 1.5 X Alkali ratio 

FIGURE 2.2.1: Alkali ratio- input parameter 




Legend 
MSG - Membership grade 
FS- First stage, SS-Second 
stage, TS- Third Stage 
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25 Y Kiln load material in tons 



FIGURE 2.2.2: Kiln load material in tons-output parameter 
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40 Z Percentage of CKD 



Legend 
MSG - Membership grade 
VL- very less, L-low, M-medium, 
H-high, VH- very high 



FIGURE 2.2.3: Percentage of CKD - output parameter 



For the terms of alkali ratio, kiln load material in tons and 
percentage of CKD we give the following membership 
functions: 



^( X ) a 



ii L (X) = (l-X)/0.5 0.5<X<1 
(X-0.5)/0.5 0.5<X<1 
(1.5-X)/0.5 1<X<1.5 

li h (X) = (X-1)/0.5 1<X<1.5 



M*) : 



(2.2.1) 



H(Y)- 



ratio in tons 



ii FS (Y) = (15-Y)/10 5<Y<15 

(Y-5)/10 5<Y<15 

(25-Y)/10 15<Y<25 

li ts (Y) = (Y-15)/10 15<Y<25 



Mss(Y) = 



(2.2.2) 
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n(z) 



percentage 
of CKD 



M z ) = 

Mm( Z )= 

M z ) = 



H VL (Z) = (10-Z)/10 0<Z<10 

Z/10 < Z < 10 

(20-Z)/10 10<Z<20 

(Z-10)/10 10<Z<20 

(30-Z)/10 20<Z<30 

(Z - 20)/10 20 < Z < 30 

(40 - Z)/10 30 < Z < 40 

H VH (Z) = (Z-30)/10 30<Z<40 



(2.2.3) 



By applying the "if ... and ... then" rules [refer 11] to the three- 
membership functions |~i(X), |~i(Y) and |a(Z) we get the 
following table of rules. 

The rules given in Table 2.2.1 read as follows : 



Table 2.2.1 



x\ 


FS 


SS 


TS 


L 


VL 


M 


H 


M 


L 


M 


H 


H 


M 


H 


VH 



For example : 

If alkali ratio is L and kiln load material in tons is FS then 
percentage of CKD is VL. If alkali ratio is H and kiln load 
material in tons is TS then percentage of CKD is VH; and so on. 
Rules of evaluation using the membership functions defined 
by the equation (2.2.1) and (2.2.2), if alkali ratio is 1.2 and kiln 
load material is 17 tons we get the fuzzy inputs as |a M (l-2) = 0.6, 
Hh(1.2) = 0.4, nss(17) = 0.8 and Urs(17) = 0.2. Induced decision 
table for percentage of CKD is as follows. 
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Table 2.2.2 



^^^\ Y 





Uss (17) = 0.8 


H TS (17) = 0.2 














Ml-2)=0.6 





Hm(Z) 


MZ) 


Hh(1.2)=0.4 





MZ) 


Hvh(Z) 



Conflict resolutions of the four rules is as follows: 

Rule 1 : If X is M and Y is SS then Z is M 
Rule 2 : If X is M and Y is TS then Z is H 
Rule 3 : If X is H and Y is SS then Z is H 
Rule 4 : If X is H and Y is TS then Z is VH 

Now, using Table 2.2.2 we calculate the strength values of the 
four rules as 0.6, 0.2, 0.4 and 0.2. Control output for the 
percentage of CKD is given in table 2.2.3. 







Table 2.2.3 




X ^^^ 





Li ss (17) = 0.8 


Ll TS (17) = 0.2 














li m (1.2)=0.6 





min{[0.6, li m (Z)]} 


min{[0.2, Li H (Z)]} 


l 1h (1.2)=0.4 





min{[0.4, Li H (Z)]} 


min{[0.2, livh(Z)]} 



To find the aggregate(agg) of the control outputs, we obtain 
the maximum of the minimum. This is given by the following 
figure 2.2.4, that is Li agg (Z) = max {min {[0.6, li m (Z)] min {[0.4, 
li h (Z)],)], min [0.2, |a vH (Z)]}. By applying the mean of 
maximum method for defuzzification that is the intersection 
points of the line li = 0.6 with the triangular fuzzy number 
li m (Z) in equation (2.2.3) we get the crisp output to be 20%. 
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40 Z Percentage of CKD 



FIGURE 2.2.4: Aggregated output and defuzzificztion for 
the percentage of CKD 



Rules of evaluation using the membership function defined by 
the equation (1) and (2), if alkali ratio is 0.5 and kiln load 
material is 5 tons we get the fuzzy inputs as u. L (0.5) = 1, u. H (0.5) 
= 0, u. rs (5) = 1 and u. ss (5) = 0. Induced decision table for 
percentage of CKD is as follows. 





Table 22A 


\ 




X ^^\ 


Hfs(5) = 1 


Mss(5) = 





Hl(0.5) = 1 


Hvl(Z) 


M z ) 





Mm(0.5 )= 


Hl(Z) 


Hm(Z ) 


















Conflict resolutions of the four rules is as follows: 

Rule 1 : If X is L and Y is FS then Z is VL 
Rule 2 : If X is L and Y is SS then Z is M 
Rule 3 : If X is M and Y is FS then Z is L 
Rule 4 : If X is M and Y is SS then Z is M. 

Now, using Table 2.2.4 we calculate the strength values of 
the four rules as 1, 0, and 0. Control output for the percentage 
of CKD is given in Table 2.2.5. 
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Table 2.2.5 



X ^^^ 


Hfs (5) = 1 


Mss(5) = 





Hl(0.5) = 1 


min {[1, Hvl(Z)]} 


min {[0, Hm(Z)]} 





Hh(0.5) = 


min {[0, n L (Z)]} 


min {[0, |i M ( z )]} 


















To find the aggregate of the control outputs, we obtain the 
maximum of the minimum. This is given by the following 
figure 2.2.5 that is u. agg (Z) = {min {1, u. V l(Z)]}, min{[0, 
Hm(Z)]}, min {[0, u. L (Z)]}. By applying the mean of maximum 
method for defuzzification that is the intersection points of the 
line (j. =1 with the triangular fuzzy number u. V l(Z) in equation 
(3) and get the crisp output to be 0%. 




40 Z Percentage of CKD 



FIGURE 2.2.5: Aggregated output and defuzzificztion for 
the percentage of CKD 

Rules of evaluation using the membership function defined 
by the equations (1) and (2), if alkali ratio is 1 and kiln load 
material is 15 tons we get the fuzzy inputs as u.l(1) = 0, u. H (1) = 
and Ml) = 1, Hfs(15) = 0, li ss (15) = 1, li T s(15) = 0, Induced 
decision table for percentage of CKD is as follows. 
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Table 2.2.6 



\ Y 
X ^\ 


Hfs(15) = 


Hss(15) = 1 


Hts(15) = 


Ml(1) = 


Mvl(Z) 


Hm(Z) 


Mh(Z) 


Ml) = i 


Hl(Z) 


Hm(Z) 


Mh(Z) 


Mi) = o 


Hm(Z) 


MZ) 


Mvh(Z ) 



Conflict resolutions of the nine rules is as follows : 



Rulel 


IfXis 


Rule 2 


IfXis 


Rule 3 


IfXis 


Rule 4 


IfXis 


Rule 5 


IfXis 


Rule 6 


IfXis 


Rule 7 


IfXis 


Rule 8 


IfXis 


Rule 9 


IfXis 



L and Y is FS then Z is VL 
L and Y is SS then Z is M 
L and Y is TS then Z is H 
M and Y is FS then Z is L. 
MandYisSSthenZisM. 
M and Y is TS then Z is H. 
H and Y is FS then Z is L. 
H and Y is SS then Z is M. 
H and Y is TS then Z is H. 



Now, using Table 2.2.6 we calculate the strength values of the 
nine rules as 0, 0, 0, 0, 1, 0, 0, 0, 0. Control output for the 
percentage of CKD is given in Table 2.2.7. 





Ta 


ble 2.2.7 






Hfs (15) = 


Hss(15) = 1 


urs(15) = 1 


Hl(1)=0 


min{[0,Li V L(Z)]} 


min{[0,Li M (Z)]} 


min{[0,Li H (Z)]} 


Hm(1)=1 


min{[0,Li L (Z)]} 


min{[0,Li M (Z)]} 


min{[0,Li H (Z)]} 


Hh(1)=0 


min{[0,LiM(Z)]} 


min{[0,Li H (Z)]} 


min{[0,Li H (Z)]} 



To find the aggregate of the control outputs, we obtain the 
maximum of the minimum. This is given by the following 
figure 2.2.6, that is |a agg (Z) = max {min {0, li V l(Z)]}, min{[0, 
MZ)]}, min {[0, Li L (Z)]}, {min {1, ii H (Z)]}, min{[0, li V h(Z)]}. 
By applying the mean of maximum method for defuzzification 
that is the intersection points of the line |a =1 with the triangular 
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fuzzy number u. V l(Z) in equation (2.2.3) and get the crisp output 
to 20%. 




40 Z Percentage of CKD 



FIGURE 2.2.6: Aggregated output and defuzzification for 
the percentage of CKD 



Rules of evaluation using the membership function defined by 
the equations (2.2.1) and (2.2.2), if alkali ratio is 1.5 and kiln 
load material is 25 tons we get the fuzzy inputs as u. M (1-5) = 0, 
u-h (1.5) = 1, u.ss (25) = and u, TS (25) = 1. Induced decision 
table for percentage of CKD is as follows. 

Table 2.2.8 



^^\ Y 





Mss(25) = 


Mts(25) = 1 














Mi-5) = o 





Hm(Z) 


MZ) 


Mi-5) = 1 





Hh(Z) 


Mvh(Z ) 



Conflict resolutions of the four rules is as follows : 

Rule 1 : If X is M and Y is SS then Z is M 
Rule 2 : If X is M and Y is TS then Z is H 
Rule 3 : If X is H and Y is SS then Z is H 
Rule 4 : If X is H and Y is TS then Z is VH. 
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Now, using Table 2.2.8 we calculate the strength values of 
the four rules as 0, 0, and 1 Control output for the percentage 
of CKD is given in Table 2.2.9. 

Table 2.2.9 



X ^^\ 





Hss(25) = 


Mts(25) = 1 














Mi-5) = o 





min [0, n M (Z)] 


min [0, Hh(Z)] 


Mh(1-5) = 1 





min [0, |i H (Z)] 


min [1, livh(Z)] 



To find the aggregate of the control outputs, we obtain the 
maximum of the minimum. This is given by the following 
figure 2.2.7, that is u. agg (Z) = max {min {0, u, M (Z)]} , min {[0, 
Mh(Z)]}, min{[l, u. V h (Z)]}. By applying the mean of maximum 
method for defuzzification that is the intersection points of the 
line (a =1 with the triangular fuzzy number u.vh(Z) in equation 
2.2.3 and get the crisp output to 40%. 




40 Z Percentage of CKD 



FIGURE 2.2.7: Aggregated output and defuzzification for 
the percentage of CKD 

From our study we suggest in the online process to reduce 
(or) minimize the amount of CKD in the industry one should 
change the condition of fuel burning system and other system in 
kiln from time to time depending on the percentage of CKD in 
tons given above. 
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2.3 Determination of gas volume setpoint and 
temperature set point for CKD processing 

The total CKD dust carried out from the kiln is again returned to 
the kiln as a feed (Recycled CKD). After recycled process, we 
get some amount of remaining CKD from kiln, which is 
disposed in the environment(as a waste polluting the 
environment). Most of the cement factory uses electrostatic 
precipitator(ESP) method for recycling process of CKD, as it 
operates by gas volume and temperature. In ESP, we mainly 
concentrate on gas volume in m 3 /minute and temperature degree 
in celsius. The range of gas volume is varying from 11865 to 
15174 m 3 /minute and temperature is varying from 350°C to 
450°C. When in the recycle; the clinker is got from the 
reproduced dust to clinker by pre heater in dust collector(ESP). 
Generally an industry to minimize the net CKD dust upto 20% 
by reprocessing method randomly chooses the gas volume and 
temperature from the range of gas volume (11865 to 15174 
m 3 /minute) and temperature (350°C to 450°C) respectively. 
Since the gas volume and temperature are main concerns on 
ESP, the reprocessing directly depends on gas volume and 
temperature. The randomly choosing of the gas volume set point 
and temperature set point from the ranges of gas volume and 
temperature is uncertain and does not usually give the desired 
outcomes so, this gas volume and temperature affect the CKD 
reprocessing largely. In order to over come these problems we 
use fuzzy control to find the set point of gas volume and 
temperature in ESP, which is described in the following. The 
ranges of gas volume, temperature and percentage of net CKD 
are measured from the past happening data in ESP on a scale, 
are 11865 to 15174 m 3 /minute, 350°C to 450°C and to 20% 
respectively. Temperature = {350°C [low (L)], 400°C [medium 
(M)] and 450°C [high (H)]}. Gas volume = 11865 to 
15174m 3 /min [first stage (FS)], 13020 m 3 /min [second stage 
(SS)], and 15174 m 3 /min [third stage (TS)}. Percentage of net 
CKD = {0[very less (VL)], 5 [less (L), 10[medium (M)], 15 
[high (H)] and 20 [very high (VH)}. The terms of these 
parameters are presented in figures 2.3.1 or 2.3.3. 
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350 400 450 X Temperature 

FIGURE 2.3.1: Temperature - input parameter 




11865 13020 15174 Y Gas Volume 
FIGURE 2.3.2: Gas Volume- input parameter 




5 10 15 20 Z Percentage of CKD 

FIGURE 2.3.3: Percentage of net CKD-output parameter 
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For the terms of temperature, gas volume and percentage of net 
CKD we give the following membership functions : 



n(x) 



temperature 



H L (X) = (400-X)/50 350<X<400 

f(X-350)/50 350<X<400 
H M (X)= '' (2.3.1) 

[(450-X)/50 400<X<450 

H H (X) = (X-400)/50 400<X<450 



M(Y) B 



n(z) 



Mss(Y) : 



percentage 
of net CKD 



li fs (Y) = (13020 - Y)/1155 11865 < Y < 13020 
(Y-11865)/1155 11865<Y<13020 
(15174 -Y)/2154 13020<Y<15174 
li ts (Y) = (Y -13020)/2154 13020<Y<15174 

(2.3.2) 
li vl (Z) = (5-Z)/5 0<Z<5 

Z/5 < Z < 5 

(10-Z)/5 5<Z<10 

(Z-5)/5 5<Z<10 

(15-Z)/5 10<Z<15 

(Z-10)/5 10<Z<15 

(20-Z)/5 15<Z<20 

li vh (Z) = (Z-15)/5 15<Z<20 



Ml(Z) = 
MmC 2 ) 
Mh( z ) 



(2.3.3) 



By applying the if . . . and . . . then rules to the three-membership 
function |~i(X), |a(Y) and |~i(Z), we get the following table of 
rules. The rules given in Table 2.3.1 read as follows: 

Table 2.3.1 



\y 
x^\ 


FS 


SS 


TS 


L 


VL 


M 


H 


M 


L 


M 


H 


H 


M 


H 


VH 
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For example: 

If temperature is L and gas volume is SS then percentage of 
net CKD is M. 

If temperature is M and gas volume is TS then percentage 
of net CKD is H; and so on. 

Rules of evaluation using the membership functions defined 
by the equation (2.3.1) and (2.3.2), if temperature is 430°C and 
gas volume is 13080 m 3 /min we get the fuzzy inputs as |a M (430) 
= 0.4, li h (430) = 0.6, li ss (13080) = 0.97 and li T s(3080) = 0.02. 
Induced decision table for percentage of net CKD is as follows. 

Table 2.3.2 



X ^^^^ 





li ss (13080)= 0.97 


Li T s(13080)= 0.02 














li m (430)=0.04 





Hm(Z) 


MZ) 


li h (430)=0.06 





Hh(Z) 


Mvh(Z ) 



Conflict resolutions of the four rules is as follows : 

Rule 1 : If X is M and Y is SS then Z is M 
Rule 2 : If X is M and Y is TS then Z is H 
Rule 3 : If X is H and Y is SS then Z is H 
Rule 4 : If X is H and Y is TS then Z is VH. 

Now, using Table 2.3.2 we calculate the strength values of 
the four rules as 0.4, 0.02, 0.06 and 0.02. Control output for the 
percentage of net CKD is given in Table 2.3.3. 

Table 2.3.3 








Liss (13080) = 0.97 


li ts (13080) = 0.02 














li m (430)=0.4 





min {[0.4, MZ)]} 


min {[0.02, ii H (Z)]} 


li h (430)=0.6 





min {[0.6, n H (Z)]} 


min{[0.02, li V h(Z)]} 
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To find the aggregate of the control outputs, we obtain the 
maximum of the minimum. This is given by the following 
figure 2.3.4, that is u. agg (Z) = max{min [0.4, u. M (Z)]}, min{[0.6, 
M-m(Z)]}, min {[0.02, u. V h (Z)]}. We apply the mean of 
maximum method for defuzzification that is the intersection 
points of the line u, = 0.6 with the triangular fuzzy number 
li h (Z) in equation (2.3.3) and get the crisp output as 15 to 20%. 




20 Z Percentage of CKD 



FIGURE 2.3.4: Terms of the output and defuzzification for 
the percentage of net CKD 

Rules of evaluation using the membership functions defined by 
the equation (2.3.1) and (2.3.2), if temperature is 350°C and gas 
volume is 11865 m 3 /min we get the fuzzy inputs as u. L (350) = 1, 
u. H (350) = 0, u. FS (11865) = 1 and u. ss (11865) = 0. Induced 
decision table for percentage of net CKD is as follows. 

Table 2.3.4 



X ^^^ 


lifs (11865) = 1 


liss(H865) = 





m,(350) = 1 


u.vl(Z) 


Hm(Z) 





li m (350) = 


Hl(Z) 


Hm(Z) 


















Conflict resolutions of the four rules is as follows : 
Rule 1 : If X is L and Y is FS then Z is VL 
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Rule 2 : If X is L and Y is SS then Z is M 
Rule 3 : If X is M and Y is FS then Z is L 
Rule 4 : If X is M and Y is SS then Z is M. 

Now, using Table 2.3.4 we calculate the strength values of 
the four rules as 1, 0, and 0. Control output for the percentage 
of net CKD is given in table 2.3.5. 

Table 2.3.5 





li fs (11865) = 1 


Li ss (11865 ) = 





li l (350)=0.4 


min {[1, livl(Z)]} 


min {[0, Li M (Z)]} 





Hm(350)= 0.6 


min {[0, li l (Z)]} 


min {[0, Hm(Z)]} 


















To find the aggregate of the control outputs, we obtain the 
maximum of the minimum. 

This is given by the following figure that is u. agg (Z) = 
{min{[l, livl (Z)]}, min{[0, li m (Z)]}, min {[0, li l (Z)]}. We 
apply the mean of maximum method for defuzzification that is 
the intersection points of the line |u = 1 with the triangular 
fuzzy number u. V l(Z) in equation (2.3.3) and get the crisp output 
as to 5%. 




20 Z Percentage of CKD 



FIGURE 2.3.5: Terms of the output and defuzzification for 
the percentage of net CKD 
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Rules of evaluation using the membership functions defined 
by the equation (4) and (5), if temperature is 400°C and gas 
volume is 13020 m 3 /min we get the fuzzy inputs as li l (400) = 0, 
li m (400) = 1, li h (400) = 0, li fs (13020) = 0, li S s(13020) = 1 and 
li T s(13020) = 0. Induced decision table 2.3.6 for percentage of 
net CKD is as follows. 





Tab] 


e 2.3.6 




\^ Y 


li fs (13020)=0 


li ss (13020)=1 


Li TS (13020)=0 


Hl(400)=0 


Hvl(Z) 


Hm(Z) 


MZ) 


Hm(400)=1 


Hl(Z) 


M z ) 


M z ) 


li h (400)=0 


Hm(Z) 


MZ) 


Hvh(Z ) 



Conflict resolutions of the nine rules is as follows: 

Rule 1 : If X is L and Y is FS then Z is VL 
Rule 2 : If X is L and Y is SS then Z is M 
Rule 3 : If X is L and Y is TS then Z is H 
Rule 4 : If X is M and Y is FS then Z is L. 
Rule 5 : If X is M and Y is SS then Z is M. 
Rule 6 : If X is M and Y is TS then Z is H. 
Rule 7 : If X is H and Y is FS then Z is M. 
Rule 8 : If X is H and Y is SS then Z is H. 
Rule 9 : If X is H and Y is TS then Z is VH. 

Now, using Table 2.3.6 we calculate the strength values of 
the nine rules as 0, 0, 0, 0, 1, 0, 0, 0, 0. Control output for the 
percentage of net CKD is given in Table 2.3.7. 





Tab 


Le 2.3.7 




~\Y 


li fs (13020) = 


li ss (13020) = 1 


Li TS (13020) = 


Hl(400) =0 


min{[l,nvL(Z)]} 


min{[0,Li M (Z)]} 


min{[0,Li H (Z)]} 


li m (400)=1 


min{[0, li l (Z)]} 


min{[l,Li M (Z)]} 


min{[0,Li H (Z)]} 


li h (400)=0 


min{[0, Li M (Z)]} 


min{[0,Li H (Z)]} 


min{[0,Liv H (Z)]} 
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To find the aggregate(agg) of the control outputs, we obtain 
the maximum of the minimum. 

This is given by the following figure 2.3.6, that is |a agg (Z) = 
max {min {[0, li V l(Z)]}, min {[1, n M (Z)],)L min {[0, |a L (Z)]}, 
min {[0, |a H (Z)]}, min {[0, li V h(Z)]}. We apply the mean of 
maximum method for defuzzification that is the intersection 
points of the line li = 1 with the triangular fuzzy number lim(Z) 
in equation (2.3.3) we get the crisp output to be 10 % to 15 %. 




20 Z Percentage of net 
CKD 

FIGURE 2.3.6: Aggregated output and defuzzification for 
the percentage of net CKD 

Rules of evaluation using the membership functions defined 
by the equation (2.3.1) and (2.3.2), if temperature is 450°C and 
gas volume is 15174 m 3 /min we get the fuzzy inputs as li m (450) 
= 0, li h (450) = 1, li ss (15174) = and ii TS (15174) = 1. Induced 
decision table for percentage of net CKD is as follows. 

Table 2.3.8 



X ^^^ 





^ss (15174) = 


li ts (15174) = 1 














li m (450) = 





Hm(Z) 


MZ) 


Hh(450)= 1 





Hh(Z) 


Mvh(Z) 



Conflict resolutions of the four rules is as follows: 
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Rule 1 : If X is M and Y is SS then Z is M 
Rule 2 : If X is M and Y is TS then Z is H 
Rule 3 : If X is H and Y is SS then Z is H 
Rule 4 : If X is H and Y is TS then Z is VH. 

Now, using Table 2.3.8 we calculate the strength values of the 
four rules as 0, 0, and 1. Control output for the percentage of 
net CKD is given in Table 2.3.9. 

Table 2.3.9 








Hss(15174) = 


Uts(15174) = 














li m (450)=0 





min[0, Um( z )] 


min[0,Li H (Z)] 


Mh(400)= 1 





min[0, |i H (Z)] 


min[l,nvH(Z)] 



To find the aggregate(agg) of the control outputs, we obtain 
the maximum of the minimum. 

This is given by the following figure 2.3.7, that is u. agg (Z) = 
max{min{[0, li m (Z)]}, min{[0, li h (Z)],)], min{[l, li V h(Z)]}. We 
apply the mean of maximum method for defuzzification that is 
the intersection points of the line |u = 1 with the triangular fuzzy 
number u.m(Z) in equation (2.3.3) we get the crisp output to be 
20 %. 




20 Z Percentage of net 
CKD 

FIGURE 2.3.7: Aggregated output and defuzzification for 
the percentage of CKD 
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2.4 Finding the MIX of raw materials in proper proportion 
and minimize the waste dust using fuzzy neural network 

The study of proper proportions of material mix during the 
clinkerization process is very difficult due to inconsistency in 
the chemical and mineralogical composition and the variation of 
these characteristic affects kiln operation, fuel consumption, 
clinker quality and above all the amount of CKD vent into the 
atmosphere. Further the raw mix should maintain a fixed range 
for a specific quality of cement. The problem of satisfying this 
range involves lot of randomness and uncertainty, which in turn 
speaks about the desired quality of the clinker. Chemical and 
mineralogical composition contains Si0 2 , A1 2 3 , Fe 2 3 , CaO, 
MgO, K 2 and Na 2 0. Since all terms used to determine the 
proper proportions of material mix is very ambiguous, we felt it 
would be proper to use fuzzy theory approach to study the 
problem. We adopt fuzzy relational neural network method to 
find the correct proportion of raw mix so that the desired quality 
of the clinker is achieved. This is done by taking experts 
opinion about the proportions and then by giving fuzzy weights. 
This membership grades are varied a finite number of times till 
the error function reaches zero, which is equivalent to studying 
the set point values. The clinker of desired chemical 
composition is expected to satisfy the following modulus related 
to the chemical composition of the raw mix. 
Lime saturation factor (LSF), 



CaO x 100 

LSF = (2.4.1) 

2.8 SiQ 2 + 1.2 A1 2 3 + 0.65Fe 2 O 3 



A high LSF requires high heat consumption for clinker 
burning inside the kiln and this gives more strength to the 
cement. 
Silica Modulus (SM) 

SM = ^2 (2.4.2) 

Al 2 3 + Fe 2 3 
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A higher SM decreases the liquid phase content, which 

impairs the burnability of the clinker and reduces the cement 

setting time. 

Alumina Modulus (AM) 

ai n 
AM = ^^ (2.4.3) 

Fe 2 3 



The value of AM determines the composition of liquid phase in 
the clinker. 

Here we describe the working of the block schematic of raw 
mill processing. The raw mill grinder receives raw materials 
such as limestone, silica, iron and bauxite for the production of 
cement in separate feeders, called weigh feeders. All the raw 
materials are ground in a raw mill grinder to a powder form. A 
sample of this ground raw mix is collected at the output of the 
raw mill grinder by an auto sampler, and a sample is prepared 
after being finely ground by vibration mill and pressed by 
hydraulic press and then is analysed in the laboratory by an X- 
ray sequential spectrometer. The results of X-ray analysis, 
which are obtained through sampling and analysing the 
equipment, are fed to the computer through a data 
communication line, for the required control action. The entire 
process is illustrated in figure 2.4.1. 

Weigh 

feeders „ . 

Auto 



V 




Raw mill 


Sampler 


Storage & 
Blending 




Plant 




\ 




\? 










\ y 




























r^ 


Sequential 

X-Ray 

Spectrometer 






{J 

Sample 

PreDaration 
















Computer 
System 























FIGURE 2.4.1 Block schematic of raw mill processing steps 
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The past researchers developed a control algorithm for raw 
mixing proportion based of singular value decomposition(SVD) 
methods. The purpose of this algorithm is to calculate the 
change in raw materials in each of the weigh feeders to achieve 
the raw mixing that is LSF, SM and AM. 

Singular value decomposition(SVD) is one of the most 
basic and important tools in the analysis and solution of the 
problems in numerical linear algebra, and are finding increasing 
applications in control and digital signal processing. The 
potential of SVD technique is first exploited in the domain of 
linear algebra, where it provides a reliable determination of the 
rank of the matrix, thereby leading to accurate solutions of 
linear equations. 

Here we adopt raw mix proportion control algorithm to our 
problem. The purpose of this algorithm is to calculate the 
change in raw materials in each of the weigh feeders to achieve 
the target value of the chemical composition ratio (or) module 
of LSF, SM and AM. 

Suppose at any instant the action of the control system gives 
rise to the composition change as dLSF', dSM' and dAM' in 
response to the required composition change as dLSF, dSM and 
dAM respectively then the total mean square error at that instant 
will be 

E = (dLSF - dLSF') 2 +(dSM - dSM') 2 + (dAM - dAM') 2 (2.4.4) 

The problem now is to minimize E with respect to the change in 
the feeder content(dw;: i = 1, 2, ..., n). Differentiating equation 
(2.4.4) with respect to dw and equating to zero, we will have 
dLSF' = dLSF, dSM' = dSM and dAM' = dAM. As mentioned 
earlier, the values of LSF, SM and AM of the raw material, 
change constantly. 

Our objective is to keep the values of LSF, SM and AM of 
the raw mix at the raw mill outlet fixed by changing the quantity 
of the raw material in the weigh feeders. So the module LSF, 
SM and AM are functions of the change in the raw material in 
different feeders. This can be represented as 
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" BJ SF 

dSLF' = dLSF = V dw i (2.4.5) 

dSM' = dSM = V dw. (2.4.6) 

dAM' = dAM = V dWj (2.4.7) 



Jdw^O (2.4.8) 

i=l 

LL ; < dw; < HL; (2.4.9) 



where w ; is the mix ratio of raw material in the feeder, LL ; and 
HL ; are the lower limit and the higher limit respectively of the 
raw material change possible for the i th feeder(i = 1, 2, ..., n). 
The composition change, for example in LSF is given by 

dLSF' = LSF sp - LSF meas (2.4.10) 

Here 'sp' stands for set point that is the desired value 
and 'meas' stands for the measured value that is the value 
achieved. 

Now consider the solution of equation (2.3.5) to (2.3.8). The 
number of unknowns is the same as the number of weigh 
feeders. If there are four unknowns then there are four weigh 
feeders, we have the following set of equations with four 
unknowns. 

SLSF_, (3LSF_, dLSF , dLSF , JT „„ 
dw, H dw, H dw, h dw. = dLSF 

ow 1 ow 2 ow 3 ow 4 

(2.4.11) 
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5SM_, SSM_, SSM_, SSM_, ,„ w , 

-dw, h dw, + dw, + dw, =dSM 



dw 1 



3w, 



5w, 



dw„ 



(2.4.12) 



5AM j 5AM , 5AM , 5AM , ,,„. 

-dWj + dw, + dw, + dw, = dAM 



dw 1 



9w, 



dw, 



dw. 



(2.4.13) 

dWj + dw 2 +dw 3 + dw 4 =0 (2.4.14) 

Rearranging equations (2.4.11) to (2.4.14) in matrix form yields 
(3LSF 5LSF 5LSF SLSF~ 



dw 1 dw 2 



dw 3 dw 4 

5sm esM esM ssm 



dw, 



5w 2 5w 3 dw 4 



5AM 5AM 5AM 5AM 



dw 1 dw 2 



dw. dw. 



1 



1 



1 



1 



dw 1 




"dLS F'~ 


dw 2 




dSM" 


dw 3 




dAM' 


|_dw 4 _ 








(2.4.15) 



If the stacker reclaimed, a macline that feeds limestone of 
constant chemical composition to the weigh feeders is available, 
then LSF value will more(or) less remain constant; So in this 
case, one must give importance to achieving desired value for 
SM and AM. To cope with this situation in SVD method one 
can simply ignore equation(2.4.11). Also this method can be 
used in the event of feeder failure, or the addition of a feeder. In 
these cases, the number of feeder is simply changed and the 
corresponding equations, similar to equation (2.4.11) are added 
(or) deleted as appropriate. 

The value is the amount of change for that modulus with 
unit change in raw material mix proportion sent into the grinder. 
This can be obtained from the calculation of the composition of 
the raw materials, but in cement production process the 
composition of the raw materials fed into the mill changes 
constantly. So it is not possible to get fixed values for these 
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differential factors. Raw materials from a particular quarry have 
the composition varying over very narrow ranges for your 
purpose we have chosen a typical composition of raw material 
with its values as the average value of the material received 
from the quarry. The raw materials in each feeder consist of 
CaO, Si0 2 , A1 2 3 , and Fe 2 3 thus affecting all the three moduli 
such as LSF, SM and AM as given in equations (2.4.1, 2.4.2 and 
2.4.3) so these moduli can now be redefined as 

n 

Z Ca °i w i 

LSF = - ^ 

X[2.8(SiO 2 )-w i + 1.2(Al 2 O 3 ) i -w i + 0.65(Fe 2 O 3 )-w i 

i=l 

(2.4.16) 

X(Si 2 3 ) 1 -w 1 
SM = — ^ (2.4.17) 

X[(Al20 3 ) 1 -w 1 + (Fe 2 3 ) 1 - w i] 



AM = ^ (2.4.18) 

SRFeAX-w,] 



where n is the number of feeders. Now the differential 
coefficients of equation (2.4.11), (2.4.12) and (2.4.13) can be 
obtained by differentiating the equation (2.4.16), (2.4.17) and 
(2.4.18) with respect to w ; . 

Adaptation of fuzzy neural network to raw mix proportion 
control algorithm : 

Let P represent the coefficient of raw mixing ratio that is 

5LSF 5SM 5AM ,.,„„,„ 

, , where 1 = 1, 2, 3, 4, Q represents the 

3wj <5W; <9wj 

unknown quantities for four weigh feeders that is dwi, dw 2 , dw 3 
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and dw 4 and R represents the known values that is dLSF', dSM' 
and dAM'. Generally researchers used some other non-fuzzy 
method to estimate the unknowns dwi, dw 2 , dw 3 and dw 4 but 
since one is not always certain of solving these equations, fuzzy 
neural network model is adopted. By this method one is always 
guaranteed of a solution. 

The problem is tackled in two stages according as if a 
solution exist using fuzzy relation equations PoQ = R then all 

the quantities for four weigh feeders are determined. If PoQ = R 

does not give solution the fuzzy neural network method is 
adapted to the fuzzy relation equation as the second stage. By 
adopting fuzzy neural network method to the fuzzy relation 
equation, unknown quantities for four weigh feeders are 
determined. 

We get the matrices according to P ° Q = R. 
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(2.4.19) 



where 
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Q = 



dw 1 




"dLS F'~ 


dw 2 

dw 3 


andR = 


dSM' 
dAM' 


_dw 4 _ 








we in this problem minimize the error between the rise to the 
composition change and required composition change. The 
membership value py e [0,1] are given by experts. 

Equation (2.4.19) can be rewritten as 



Pll Pl2 Pl3 Pl4 

P21 P22 P23 P24 

P31 P32 P33 P34 

P41 P42 P43 P44 



dw : 
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It can partitioned into 4 equations. 
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and 



[P 4 



P44] 



dw 1 




"dLS F'~ 


dw 2 




dSM' 


dw 3 




dAM' 


_dw 4 _ 








If the above partitioned equation do not satisfy the condition 
max q ik < r ik (where q ik are unknown quantities for weigh 
feeders and r ik are known values that is dLSF', dSM' and dAM') 
then the system of equations has final solution. If the above 
partitioned equation satisfy this condition max q ik < r ik where q ik 
are unknown quantities for weigh feeders and r ik are known 
values that is dLSF', dSM' and dAM' then the system of 
equations has no solution. In this case fuzzy neural network 
method is adopted for fuzzy relation equation as the second 
stage. 
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dw : 
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dw 2 
dw 3 


= 


dSM' 
dAM' 
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The linear activation function f defined earlier gives the output 
yi=f[max (wy Xj)], i e N n . First calculate w u Xi, Wi 2 x 2 , Wi 3 x 3 , and 
W14X4 then find yi= f [max(wij Xj)] which gives dwi. Similarly 
calculate w 2 iXi, w 2 2X 2 , w 23 x 3 , and w 24 x 4 to find y 2 =f [max (w 2 j 
Xj)] which gives dw 2 , calculate w 3 iXi, w 32 x 2 , w 33 x 3 , and w 34 x 4 to 
find y 3 =f [max (w 3j Xj)] which gives dw 3 and calculate w 4 iX 1; 
w 42 x 2 , w 43 x 3 , and w 44 x 4 to find y 4 =f [max (w 4 j Xj)] which gives 
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dw 4 . Then we can find out using equation (2.4.4) that is whether 
the error function reaches or not, suppose the error function 
does not reach 0, we change the weights that is the membership 
grades till the error reaches zero, which is explained in 
figure(2.4.8). Thus to achieve the value of error function to be 
zero, we give different membership grades to the weigh feeders 
(finite number of times) and make the value of required 
composition change to be equal to the raise in composition 
change. 



Xi 


Wn 


W 2 1 


W31 


W41 


X2 




Wl2 




W 2 2 




W32 




W42 


X3 






Wl3 






W 2 3 






W33 






W43 


x 4 








Wl4 








W 24 








W34 











ONi 



ON, 



ON 3 



W 44 



ON 4 



Yj Y 2 Y 3 Y 4 

Figure 2.4.2: The feed forward neural network 



2.5 Conclusions 

The fuzzy control method described and defined in this chapter 
has the following problem: 

1) Monitoring and control of the system 

2) CKD reprocessing 

Monitoring and control : 

We have analyzed the alkali ratio (0.5 to 1.5) in kiln load 
material. The alkali ratio and the kiln load material in tons are 
considered as the two input parameter of fuzzy control. The 
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output, is the percentage of CKD. The estimated results are as 
follows : 

1) When the alkali ratio is 0.5 in the 5 tons of kiln load 
material, the estimated percentage of CKD is 0. 

2) When the alkali ratio is 1 in the 15 tons of kiln load 
material, the estimated percentage of CKD is 20. 

3) When the alkali ratio is 1.2 in the 17 tons of kiln load 
material, the estimated percentage of CKD is 20. 

4) When the alkali ratio is 1.5 in the 25 tons of kiln load 
material, the estimated percentage of CKD is 40. 

From our study we suggest in the online process to reduce 
or minimize the amount of CKD in the industry one should 
change the condition of fuel burning system and other system in 
kiln from time to time depending on the percentage of CKD in 
tons. 

CKD Reprocessing : 

CKD reprocessing, mainly concentrates on gas volume and 
temperature set point using fuzzy control. The fuzzy control 
method suggests the following results to minimize the CKD in 
reprocessing. 

1) The suggested gas volume is 11865 m 3 /min and temperature 
set point 350°C for reprocessing of CKD. At the time of 
reprocessing with suggested gas volume and temperature set 
point, the percentage of net CKD occurs from to 5. 

2) The suggested gas volume is 13020 m 3 /min and temperature 
set point 400°C for reprocessing of CKD. At the time of 
reprocessing with suggested gas volume and temperature set 
point, the percentage of net CKD occurs from 10 to 15. 

3) The suggested gas volume is 13080 m 3 /min and temperature 
set point 430°C for reprocessing of CKD. At the time of 
reprocessing with suggested gas volume and temperature set 
point, the percentage of net CKD is 20. 

4) The suggested gas volume is 15174 m 3 /min and temperature 
set point 450°C for reprocessing of CKD. At the time of 
reprocessing with suggested gas volume and temperature set 
point, the percentage of net CKD is 20. 
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Raw material mix using fuzzy network : 

The fuzzy neural network method defined has the following 
merits. The solution exists for all unknown weigh feeders and 
made the error between raise to the composition change and 
required composition change of the raw material close to zero. 
That means the change in raw materials in each of weigh 
feeders dwi, dw 2 , dw 3 , dw 4 is achieved by the membership 
grade. This is very important one in cement industries to 
produce a desired quality of clinker. 

The merits of fuzzy control method : 

1) In earlier method the cement industry estimated the 
percentage of CKD approximately upto 40%. The industry 
did not know how much percentage of CKD occurs in each 
process. Using fuzzy control method, the estimated 
percentage of CKD in each process. By using this, the 
industry can change some internal condition of kiln and 
minimize the CKD in the online process. 

2) The earlier methods adopted by the cement industries, 
choose the temperature set point and gas volume randomly 
for electrostatic precipitator to minimize the net CKD in 
reprocessing. But the random choice did not in general give 
the desired out comes. Using fuzzy control method, gives 
exact temperature set point and set point of gas volume 
from the range of temperature set point and gas volume for 
electrostatic precipitator. By using this temperature set point 
and gas volume set point, the industry will get the desired 
outcomes. 

The merits of the fuzzy neural networks : 

1) Solution exists to all unknown weigh feeders. 

2) The target value of the chemical composition is achieved by 
minimizing the error between raise to the composition 
change and required composition change. 
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Chapter Three 



Determination of Temperature Set 
Points for Crude Oil 



3.1 Introduction 

Study of Temperature set point in Chemical industries happen 
to be an important feature. Here we give an illustration how 
fuzzy control method is adopted for finding precise temperature 
set point to distil different crude in an oil refinery. Oil that 
comes from the ground is called the "crude oil". By cooking, the 
crude is converted to useful oil. Here the temperature set point 
plays a vital role at the time of cooking the crude oil. Since the 
quality and quantity of the crude is dependent on the 
temperature set point, the crude oil refinery has different 
temperature, set points to distil different crudes. Here, this 
chapter tries to determine a precise temperature set point for the 
crude oil refinery to maximise the distillation of the crude and 
the quantity of the crude for long hours. This chapter six 
sections. 

This study is significant because most of the crude oil 
refineries have common type of operating systems. The analysis 
of this study is focussed on Kalundborg Refinery [Ebbesen 
(1992)]. Here we approach the problem of finding the precise 
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temperature set point for different crudes using fuzzy control 
theory. The data is taken from Kalundborg Refinery [Ebbesen 
(1992)]. 

In 1995, Friedman developed a Mass and Enthalpy balance 
method and used it to improve the quality of crudes. In 1992, 
Ebbesen studied about the crude operating in Kalundborg 
Refinery. He made some derivations from the theory of 
Friedman. Finally he gave a range of temperature set points for 
the distillation of different crudes. However at the end of his 
study he made it clear that in the case of kerosene, 90% stayed 
within 1°C of its set point of temperature, in the case of naphtha 
95% distillation stayed within 1°C of its set point of 
temperature. After two hours, the quality during crude switches 
was different indicating a lower quality. 

Here we establish the result using the data taken from 
Ebbesen(1992). This gives temperature set point for the 
distillation of kerosens, naphtha and gasoil. The range of 
temperature set points and the various percentage of distillation 
are converted into the fuzzy control theory. Here membership 
grade is assigned to each temperature set point and percentage 
of distillation and fuzzy control rule is used to each temperature 
set point and percentage of distillation. Finally center max-min 
rule is used, to find the precise temperature set point for 
kerosene, naphtha and gasoil. 



3.2 Description of Crude Oil Refineries 

Using the data available from Ebbesen(1992) of the Kalundborg 
oil refinery, we analyse the data via fuzzy rules and membership 
grades of fuzzy control theory method and find the precise 
temperature set points for different crudes to maximize the 
quality and distillation of crude for long hours. Crude oil 
refinery selects temperature set points randomly from the range 
of temperature set points for the distillation of the crude; as a 
result, the quality and quantity of the processed crude are 
maintained only for very few hours. Thus, to be more precise 
the aim to find the precise temperature set points for kerosene, 
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naphtha and gasoil using the data of Ebbesen (1992) from the 
Kalundborg oil refinery. 

Kalundborg oil refinery operates with different crudes on a 
regular basis. Here this crude oil refinery distils kerosene, 
naphtha and gasoil. A schematic diagram of various streams is 
shown in the following figure [Ebbessen(1992)]. 
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FIGURE 3.2.1: CRUDE REFINERY 



where TPA - denotes the top pump-around, MPA - denotes mid 
pump-around and the BPA-bottom pump-around respectively. 
Q-denotes the heat removed. These are mainly used for 
controlling the temperature. 

The random choice of temperature set points for different 
crudes with distillation taken from Ebbesen(1992) are described 
for kerosene, naphtha and gasoil. 
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Range of temperature set points for kerosene with the 
percentage of distillation 

The crude oil refinery gives 227° C-temperature set point for 
distillation of kerosene. The crude oil refinery selects randomly 
this 227° C temperature set point from the given range of set 



points {220° C, 221° C, 222° C, 



223° C, 224° C, 225° C, 226° C, 



227° C, 228° C, 229° C, 230° C}. This temperature set point 227° 
C gives 90% distillation and it says within 1° C of its set point of 
temperature. The temperature graph is given below. 

Graph 3.2.1: Graph depicting the 90% distillation of Kerosene 
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X axis- time in minute and Y axis set point of temperature for 
Kerosene. Set point was 227°C 



Legend 
TSP: Temperature set points 



From this Ebbesen (1992) concludes that for the set point 227° 
C the distillation of kerosene was 90%. 

Range of temperature set point of naphtha with percentage 
of distillation 

The crude oil refinery gives 160° C-temperature set point for the 
distillation of naphtha. The crude oil refinery selects randomly 
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this 160° C temperature set point from the given range of set 
points {155° C, 156° C, 157° C, 158° C, 160° C, 161° C, 162° C, 
163° C, 164° C, 165° C}. This temperature set point 160° C gives 
95% distillation and it stays within 1° C of its set point of 
temperature. 

The temperature graph is given below. 

Graph 3.2.2: Graph depicting the 95% distillation of Naphtha 
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X axis- time in minute and Y axis set point of temperature for 
naphtha. Set point was 160°C 



From this Ebbesen(1992) concludes that for the set point 160° C 
the distillation of naphtha was 95%. 

Range of temperature set point of gasoil with the percentage 

of distillation 

The temperature set point -4.5° C gives 95% distillation and it 
stays within 1° C of its set point temperature. 

The temperature graph is given below. 
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Graph 3.2.3: Graph depicting the 95% distillation of gasoil 
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X axis- time in minute and Y axis set point of temperature for 
naphtha. Set point was -4.50°C 



From this Ebbesen(1992) concludes that for the set point -4.5° C 
the distillation of gasoil was 95%. 



3.3 Determination of Temperature Set-Point of Kerosene 
Resulting in Better Distillation Using Fuzzy Control Theory 

The given possible ranges of temperature set points are {220° C, 
221° C, 222° C, 223° C, 224° C, 225° C, 226° C, 227° C, 228° C, 
229° C, 230° C} and possible percentages of distillation are 
(88%, 89%, 90%, 91%, 92%} in case of kerosene as observed 
by Ebbesen(1992). As fuzzy control theory is the tool adaptable 
only when the past performance data is available, now this 
chapter considers the given possible range of temperature set 
points and distillation as the inputs of fuzzy control theory. To 
identify the precise temperature set points from the possible 
range of temperature set points, this chapter assigns membership 
grades to each input of fuzzy control theory. Here, the fuzzy 
control theory is used to find a precise temperature set point for 
kerosene. 
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In the procedure developed here membership grade is from 
the interval [0, 1] to the input of each temperature set point and 
each percentage of distillation. After membership grades are 
assigned to each input of temperature set points, the following 
graph results representing the membership grades of 
temperatures set point. 



223°C 225°C 
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temperature 
FIGURE 3.3.1 : Membership function of temperature set point 



After membership grades are assigned to each input of 
distillation, the following graph results representing the 
membership grades of distillation. 
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D: Distillation 



FIGURE 3.3.2: Membership function of percentage of distillation 

The membership grade varies from to 1. 

For getting precise temperature set point for kerosene the 
throttle variables(The grade of membership) are qualified into 
five subsets. Here fuzzy rules are used to find the possible 
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percentage of distillation for each temperature set point and the 
Center Max-Min rule is used to find a throttle membership 
grade for the existing fuzzy rules. To get the grade of 
membership to each existing fuzzy rule throttle variables are 
qualified into five subsets as follows: 
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FIGURE 3.3.3: Throttle values 
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T: Temperature set point D: Percentage of Distillation 
FIGURE 3.3.4:Fuzzy rules for the temperature set point 223°C 
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The following are the fuzzy rules : 

Rule -1 : If T is 223° C temperature set point AND D is 87% 

THEN throttle is P3. 

Rule -2 : If T is 223° C temperature set point AND D is 88% 

THEN throttle is P2. 

Rule -3 : If T is 223° C temperature set point AND D is 89% 

THEN throttle is Z. 

Rule -4 : If T is 223° C temperature set point AND D is 90% 

THEN throttle is N2. 

We conclude the throttle value to the temperature set point for 
223° C, by the above stated rules, only rule -2, and rule-3 are 
applicable that is distillation is 88% and 89% respectively. 

Rule-2 

The throttle value to the temperature set point 223° C for 88% 
distillation is calculated using figures 4.2 and 4.3. 

Throttle = (0.41+0.38)/2 = 0.395 

The graphical representation of the membership grade of the 
temperature set point 223° C for 88 percentage of distillation is 
as follows. 

Graph 3.3.1: Graphical representation of Rule 2 
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Graph 3.3.2: The two outputs are then defuzzified 
by center max-min rule 
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Rule -3 

The throttle value to the temperature set point 223° C for 89% 
distillation from figures 4.2 and 4.3 is as follows 

Throttle=(0.47+.45)/2=0.46 
The graphical representation of the membership grade of the 
temperature set point 223° C for 88% percentage of distillation 
is as follows: 

Graph 3.3.3: Graphical representation of Rule 3 
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Graph 3.3.4: The two outputs are then defuzzified 
by center max-min rule 
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Here, the Center Max-Min rule is used to find a precise 
temperature set point. 

Using Center Max-Min rule to find precise temperature set 
point for kerosene 

Throttle(grade of membership) = m(P3) x Location(P2) + m(Z) 
+ Location(N2) = 0.427. Graphs for the other rules have not 
been given explicitly but after calculations, the values are given 
as the same procedure is adopted. 
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FIGURE 3.3.5:Fuzzy rules for the temperature set point 225°C 
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The fuzzy rule for distillation of kerosene: 

Rule -1 : If T is 225°C temperature set point AND D is 88% 

THEN throttle is P3. 

Rule -2 : If T is 225°C temperature set point AND D is 89% 

THEN throttle is P2. 

Rule - 3 : If T is 225°C temperature set point AND D is 90% 

THEN throttle is Z. 

Rule - 4 : If T is 225°C temperature set point AND D is 91% 

THEN throttle is N2. 

Rule -2 

The throttle value to the temperature set point 225°C for 89% 
distillation is calculated as follows. 

Throttle= (.47+.45)/2=0.46, 
Rule -3 

The throttle value to the temperature set point 225°C for 90% 
distillation is calculated as follows. 

Throttle= (.51+.23)/2=0.37, 

Using Center Max-Min rule to find precise temperature set 
point for kerosene 

Throttle(grade of membership) = m(P3) x Location(P2) + 

m(Z) + Location(N2) 
= 0.415, 
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FIGURE 3.3.6:Fuzzy rules for the temperature set point 227°C 
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The fuzzy rule for distillation of kerosene: 

Rule -1 : If T is 227°C temperature set point AND D is 89% 

THEN throttle is P3. 

Rule -2 : If T is 227°C temperature set point AND D is 90% 

THEN throttle is P2. 

Rule - 3 : If T is 227°C temperature set point AND D is 91% 

THEN throttle is Z. 

Rule - 4 : If T is 227°C temperature set point AND D is 92% 

THEN throttle is N2. 

Here we calculate the throttle value to the temperature set 
point for 227° C, by the above stated rules, only rule -2, and 
rule-3 are applicable that is distillation is 90% and 91% 
respectively. 

RuIe-2 

The throttle value to the temperature set point 227° C for 90% 
distillation is calculated as follows: 

Throttle = (0.51+0.23)/2 
= 0.37, 

Rule-3 

The throttle value to the temperature set point 227° C for 91% 
distillation is calculated as follows. 

Throttle = (0.17+0.59)/2 
= 0.38, 

Using Center Max-Min rule to find precise temperature set 
point for kerosene 

Throttle(grade of membership) = m(P3) x Location(P2) + 

m(Z) + Location(N2) 
= 0.375, 
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FIGURE 3.3.7:Fuzzy rules for the temperature set point 229°C 

Here we obtain the throttle value to the temperature set 
point for 229° C, by the above stated rules, only rule -2, and 
rule-3 are applicable that is distillation is 91% and 92% 
respectively. 

Rule-2 

The throttle value to the temperature set point 229° C for 91% 
distillation is calculated as follows. 

Throttle = (0. 17 + 0.59)/2 = 0.38, 

Rule-3 

The throttle value to the temperature set point 229° C for 92% 
distillation is calculated as follows. 

Throttle = (0.4 + 0.4)/2 = 0.4, 

Using Center Max-Min rule to find precise temperature set 
point for kerosene 

Throttle(grade of membership) = m(P3) x Location(P2) + 

m(Z) + Location(N2) 
= 0.420, 

We have taken the range of temperature set points for 
distillation of kerosene from the crude oil refinery [Ebbesen 
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(1992)] to find a precise temperature set point. This data is 
analysed with rules of fuzzy control theory. The fuzzy rules 
expressed in terms of degree of membership grade to each 
temperature set point. Finally the ultimate membership grade 
was obtained using Center Max-Min rule for the distillation of 
kerosene. 

It has been observed that the highest membership grade 
using Center Max-Min rule was given to the temperature set 
point 223°C. 



3.4 Determination of Temperature Set Point of Naphtha 
Resulting in Better Distillation using Fuzzy Control Theory 

The given possible ranges of temperature set points are {155°C, 
156°C, 157°C, 159°C, 160°C, 161°C, 162°C, 163°C, 164°C, 
165°C} and possible distillation are {93%, 94%, 95%, 96%, 
97%} in the case of naphtha as observed by Ebbesen(1992). 
Using this data as inputs of fuzzy control theory, we identify the 
precise temperature set points from possible range of 
temperature set points. Now membership grade is assigned to 
the input of each temperature set point and each percentage of 
distillation. The following graph is represents the membership 
grades of temperature set point. 
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FIGURE 3.4.1 : Membership grade of temperature set points 



After assigning membership grades in the interval [0,1] to each 
input of the percentage of distillation the following graph is 
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obtained representing the membership grades of percentage of 
distillation. 
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FIGURE 3.4.2: Membership function of percentage of distillation 

The membership grade varies from to 1. 

For getting precise temperature set point for naphtha the 
throttle variables(the grade of membership) quantified into five 
subsets. Here fuzzy rules are used to find the possible 
percentage of distillation for each temperature set point and the 
Center Max-Min rule is used to find a throttle membership 
grade for the existing fuzzy rules. To get the grade of 
membership to each existing fuzzy rule, throttle variables are 
quantified into five subsets as follows: 
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P3: very Big Positive 

FIGURE 3.4.3: Throttle values 
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IF T is 156°CANDDis 
95% THEN throttle is N2 























FIGURE 3.4.4: Fuzzy rules for the temperature set point 156°C 

The fuzzy rule for distillation of naphtha : 

Rule -1 : If T is 156°C temperature set point AND D is 92% 

THEN throttle is P3. 

Rule -2 : If T is 156°C temperature set point AND D is 93% 

THEN throttle is P2. 

Rule - 3 : If T is 156°C temperature set point AND D is 94% 

THEN throttle is Z. 

Rule - 4 : If T is 156°C temperature set point AND D is 95% 

THEN throttle is N2. 

We calculate the throttle value to the temperature set point 
for 156°C. In the above stated rules, only rule -2, and rule-3 are 
applicable that is only we get 93% and 94% of distillation 
respectively. 

RuIe-2 

The throttle value to the temperature set point 156°C for 93% 
distillation is calculated as follows. 

Throttle = (0.34 + 0.57)12 = 0.455, 
The graphical representation of the membership grade of the 
temperature set point 156°C for 93 percentage of distillation is 
as follows. 
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Graph 3.4.1: Graphical representation of Rule 2 
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Graph 3.4.2: The two outputs are then defuzzified 
by center max-min rule 
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Legend 

MSG: Membership grade 
■ ■ : 0.34 Membership grade 
• • : 0.57 Membership grade 



Rule-3 

The throttle value to the temperature set point 156°C for 94% 
distillation is calculated as follows. 

Throttle = (0.59+0.42)72 = 0.505, 
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The graphical representation of the membership grade of the 
temperature set point 156°C for 94 percentage of distillation is 
as follows. 

Graph 3.4.3: Graphical representation of Rule 3 
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Graph 3.4.4: The two outputs are then defuzzified 
by center max-min rule 
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Here, the Center Max-Min rule is used to find a precise 
temperature set point by grade f membership(throttle)value. 
Here, the Center Max-Min rule is used to find a precise 
temperature set point. 
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Using Center Max-Min rule to find precise temperature set 
point for naphtha 



Throttle(grade of membership) 



= m(P3) x Location(P2) + 

m(Z) + Location(N2) 
= 0.730. 
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95% THEN throttle is Z 
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IFTis158°CANDDis 
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FIGURE 3.4.5: Fuzzy rules for the temperature set point 158°C 

The fuzzy rule for distillation of naphtha : 

Rule -1 : If T is 158°C temperature set point AND D is 93% 

THEN throttle is P3. 

Rule -2 : If T is 158°C temperature set point AND D is 94% 

THEN throttle is P2. 

Rule - 3 : If T is 158°C temperature set point AND D is 95% 

THEN throttle is Z. 

Rule - 4 : If T is 158°C temperature set point AND D is 96% 

THEN throttle is N2. 

We calculate the throttle value to the temperature set point 
for 158°C. In the above stated rules, only rule -2, and rule-3 are 
applicable that is distillation is 94% and 95% respectively. 

Rule-2 

The throttle value to the temperature set point 158°C for 94% is 
calculated as follows: 
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Throttle = (0.59+0.42)/2 = 0.505, 

Rule-3 

The throttle value to the temperature set point 158°C for 95% 
distillation is calculated as follows. 

Throttle = (0.5+0.35)/2 = 0.425, 

Using Center Max-Min rule to find precise temperature set 
point for Naphtha 

Throttle(grade of membership) = m(P3) x Location(P2) + 

m(Z) + Location(N2) 
= 0.717, 
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IFTis160°CANDDis 
94% THEN throttle is P3 




156°C 
158°C 
160°C 
162°C 
164°C 


















IFTis160°CANDDis 
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FIGURE 3.4.6: Fuzzy rules for the temperature set point 160°C 

The fuzzy rule for distillation of naphtha : 

Rule -1 : If T is 160°C temperature set point AND D is 94% 

THEN throttle is P3. 

Rule -2 : If T is 160°C temperature set point AND D is 95% 

THEN throttle is P2. 

Rule - 3 : If T is 160°C temperature set point AND D is 96% 

THEN throttle is Z. 

Rule - 4 : If T is 160°C temperature set point AND D is 97% 

THEN throttle is N2. 
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We obtain the throttle value to the temperature set point for 
160°C. In the above stated rules, only rule -2, and rule-3 are 
applicable that is distillation is 95% and 96% respectively. 

Rule-2 

The throttle value to the temperature set point 160°C for 95% 
distillation is calculated as follows. 

Throttle = (0.5+0.35)/2 = 0.425, 
Rule-3 

The throttle value to the temperature set point 160°C for 96% 
distillation is as follows. 

Throttle = (0.59+0.5)/2 = 0.55, 

Using Center Max-Min rule to find precise temperature set 
point for naphtha 

Throttle(grade of membership) = m(P3) x Location(P2) + 

m(Z) + Location(N2) 
= 0.7625, 
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IF T is 162°CANDDis 
98% THEN throttle is N2 























FIGURE 3.4.9: Fuzzy rules for the temperature set point 162°C 



We calculate the throttle value to the temperature set point 
for 162° C. In the above stated rule only rule -2, and rule-3 are 
applicable that is distillation is 96% and 97% respectively. 
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Rule-2 

The throttle value to the temperature set point 162° C for 96% 
distillation is calculated as follows. 

Throttle = (0.5+0.6)/2 = 0.55 
Rule -3 

The throttle value to the temperature set point 162° C for 97% 
distillation is calculated as follows 

Throttle = (0.45+.55)/2 = 0.48. 

Using Center Max-Min rule to find precise temperature set 
point for kerosene 

Throttle(grade of membership) = m(P3) x Location(P2) + 

m(Z)+Location(N2) 
Throttle = 0.79, 

We have taken a range of temperature set points for the 
distillation of naphtha from the crude oil refinery to find a 
precise temperature set point. This data analyzed with rules of 
fuzzy control theory. The fuzzy rules are expressed in terms of 
the degree of membership grade to each temperature set point. 
Finally the ultimate membership grade was obtained using the 
centre max-min rule for the distillation of naphtha. 

We observe that the highest membership grade using Center 
Max-Min rule result in maximum distillation of naphtha for the 
temperature set point 162°C. 



3.5 Determination of Temperature Set-Point of Gasoil 
Resulting in Better Distillation Using Fuzzy Control Theory 

The given possible ranges of temperature set points are 
{-5.50°C, -5.00°C, -4.50° C, -4.00° C, -3.50°C} and possible 
percentages of distillation are (93%, 94%, 95%, 96%, 97%} in 
the case of gasoil as observed by Ebbesen(1992). Using this 
data as input we use fuzzy control theory to find a precise 
temperature set point for gasoil. 

Now membership grade is assigned to the input of each 
temperature set point and each distillation. After membership 
grades are assigned to each input of temperature set points, the 
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following graph results representing the membership grades of 
temperatures set point. 



-5.50°C -5.00°C 



T T 2 



-4.50°C 



-4.00°C -3.50°C 




T 5 T 6 T 7 



Ts T9 T10 

temperature 
FIGURE 3.5.1 : Membership grade of temperature set points 



After membership grades are assigned to each input of 
temperature set points, the following graph results representing 
the membership grades of temperatures set point. 
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FIGURE 3.5.2: Membership grades of percentage of distillation 

The membership grade varies from to 1. 
For getting the precise temperature set point for gasoil the 
throttle variables(the grade of membership) are quantified into 
five subsets. Here fuzzy rules are used to find the possible 
percentage of distillation for each temperature set point and the 
Centre Max-Min rule is used to find a throttle membership 
grade for existing fuzzy rules. To get the grade of membership 
to each existing fuzzy rule throttle variables are quantified into 
five subsets as follows. 
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FIGURE 3.5.3: Throttle values 
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IFTis-5.50°CANDDis 
95% THEN throttle is N2 























FIGURE 3.5.4: Fuzzy rules for the temperature set point -5.50°C 

The fuzzy rule for distillation of gasoil : 

Rule -1 : If T is -5.50 degree Celsius temperature AND D is 92% 

THEN throttle is P3. 

Rule -2 : If T is -5.50 degree Celsius temperature AND D is 

93% THEN throttle is P2. 

Rule - 3 : If T is -5.50 degree Celsius temperature AND D is 

94% THEN throttle is Z. 

Rule - 4 : If T is -5.50 degree Celsius temperature AND D is 

95% THEN throttle is N2. 



71 



We find the throttle value to the temperature set point - 
5.50°C. In the above stated rules, only rule -2, and rule-3 are 
applicable that is only we get 93% and 94% of distillation 
respectively. 

Rule-2 

The throttle value to the temperature set point -5.50°C for 93% 
distillation is calculated as follows. 

Throttle = (0.3+0.4)/2=0.35, 
The graphical representation of the membership grade of the 
temperature set point -5.50°C for 93 percentage of distillation is 
as follows. 

Graph 3.5.1: Graphical representation of Rule 2 
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Graph 3.5.2: The two outputs are then defuzzified 
by center max-min rule 
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■ ■ : 0.3 Membership grade 
• • : 0.4 Membership grade 
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Rule-3 

The throttle value to the temperature set point -5.50°C for 94% 
distillation is as follows. 

Throttle = (0.2 + 0.6)/2 = 0.40, 
The graphical representation of the membership grade of the 
temperature set point -5.50°C for 94 percentage of distillation is 
as follows. 

Graph 3.5.3: Graphical representation of Rule 3 
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Graph 3.5.4: The two outputs are then defuzzified 
by center max-min rule 
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Here, the Centre Max-Min rule is used to find a precise 
temperature set point by grade of membership(throttle)value. 
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Using Center Max-Min rule to find precise temperature set 
point for gasoil 



irottle(grade of membership) = i 

j 
= 
= 


ti(P3) x Location(P2) 
ti(Z) + Location(N2) 
3.35 x 0.5 + 0.4 x 0.5 
3.375 


+ 
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IFTis-5.00°CANDDis 
96% THEN throttle is N2 



















FIGURE 3.5.5: Fuzzy rules for the temperature set point -5.00°C 

The fuzzy rule for the distillation of gasoil : 

Rule -1: If T is -5.00 degree Celsius temperature AND D is 93% 

THEN throttle is P3. 

Rule -2: If T is -5.00 degree Celsius temperature AND D is 94% 

THEN throttle is P2. 

Rule -3: If T is -5.00 degree Celsius temperature AND D is 95% 

THEN throttle is Z. 

Rule -4: If T is -5.00 degree Celsius temperature AND D is 96% 

THEN throttle is N2. 

We calculate the throttle value to the temperature set point 
for -5.00° C. In the above stated rules, only rule -2, and rule-3 
are applicable that is distillation is 94% and 95% of distillation 
respectively. 

Rule-2 

The throttle value to the temperature set point -5.00° C for 94% 
distillation is as follows. 

Throttle = (0.6 + 0.5)/2 = 0.55 
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Rule-3 

The throttle value to the temperature set point -5.00°C for 95% 
distillation is as follows. 

Throttle = (0.35 + 0.4)/2 = 0.375, 

The graphical representation of the membership grade of the 
temperature set point -5.50°C for 94 percentage of distillation is 
as follows. 

Using Center Max-Min rule to find precise temperature set 
point for gasoil 

Throttle(grade of membership) = m(P3) x Location(P2) + 

m(Z) + Location(N2) 
Throttle = 0.55 x 1 + 0.375 x 0.5 / 0.55 + 0.375 
= 0.789. 
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IFTis-4.50°CANDDis 
97% THEN throttle is N2 























FIGURE 3.5.6: Fuzzy rules for the temperature set point -4.50°C 

The fuzzy rule for the distillation of gasoil : 

Rule -1 : If T is -4.50 degree Celsius temperature AND D is 94% 

THEN throttle is P3. 

Rule -2 : If T is -4.50 degree Celsius temperature AND D is 95% 

THEN throttle is P2. 

Rule -3 : If T is -4.50 degree Celsius temperature AND D is 96% 

THEN throttle is Z. 
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Rule -4 : If T is -4.50 degree Celsius temperature AND D is 97% 
THEN throttle is N2. 

We calculate the throttle value to the temperature set point 
for -4.50° C. In the above stated rules, only rule -2, and rule-3 
are applicable that is distillation is 95% and 96% of distillation 
respectively. 

RuIe-2 

The throttle value to the temperature set point -4.50° C for 95% 
distillation is calculated as follows. 

Throttle = (0.35 + 0.4)/2 = 0.375 

Rule-3 

The throttle value to the temperature set point -4.50°C for 96% 
distillation is calculated as follows. 

Throttle = (0.3 + 0.25)/2 = 0.277, 

Using Center Max-Min rule to find precise temperature set 
point for gasoil 

Throttle(grade of membership) = m(P3) x Location(P2) + 

m(Z)+Location(N2) 
Throttle = 0.375 x 0.5 + 0.277 x 0.5 / 0.375 + 0.277 
= 0.5 
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FIGURE 3.5.7: Fuzzy rules for the temperature set point -4.00°C 
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The fuzzy rule for the distillation of gasoil : 

Rule -1 : If T is -4.00 degree Celsius temperature AND D is 95% 

THEN throttle is P3. 

Rule -2 : If T is -4.00 degree Celsius temperature AND D is 96% 

THEN throttle is P2. 

Rule -3 : If T is -4.00 degree Celsius temperature AND D is 97% 

THEN throttle is Z. 

Rule -4 : If T is -4.00 degree Celsius temperature AND D is 98% 

THEN throttle is N2. 

Consider the temperature is -4.00 degree Celsius and the 
distillation of gasoil being 96% and 97%. Here rule-1 and rule-4 
are not applicable 

RuIe-2 

The throttle value to the temperature set point -4.00° C for 96% 
distillation is calculated as follows. 

Throttle = (0.25 + 0.3)/2 = 0.275 

Rule-3 

The throttle value to the temperature set point -4.00°C for 97% 
distillation is calculated as follows. 

Throttle = (0.4+0.4)/2 = 0.4, 

Using Center Max-Min rule to find precise temperature set 
point for gasoil 

Throttle(grade of membership) = m(P3) x Location(P2) + 

m(Z)+Location(N2) 
Throttle = 0.275 x 0.5 + 0.4 x 0.5 / 0.275 + 0.4 
= 0.501. 

We have taken a range of temperature set points for 
distillation of gasoil from the crude oil refinery [Ebbesen 
(1992)] to find a precise temperature set point. This data is 
analysed with rules of fuzzy control theory. The fuzzy rules are 
expressed in terms of degree of membership grade to each 
temperature set point. Finally the ultimate membership grade is 
obtained using Center Max-Min rule for the distillation of 
gasoil. 
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The authors have observed that the highest membership 
grade for -5.00°C using center max-min rule results in 
maximum distillation of gasoil and gives better quality. 



3.6 Conclusions 

Finding of precise set point temperatures for the distillation of 
kerosene, naphtha and gasoil have always remained to be 
uncertain in a crude oil refinery. Fuzzy control theory is able to 
predict the precise set point of temperature for kerosene, 
naphtha and gasoil, which guarantees the maximum percentage 
of distillation, and also the quality for long hours. By this 
method the random choice of temperature set point from the 
range of temperature set points, which affects the quality and 
quantity of crude is completely over come. 
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Chapter Four 



Study of Flow Rates 
in Chemical Plants 



This chapter has 3 sections. Use of FRE to estimate flow rates 
in chemical plants forms the section one of this chapter. In 
section two fuzzy neural networks are used to estimate velocity 
of flow distribution in a pipe network. The final section 
estimates the three-stage counter current extraction unit again 
using fuzzy neural networks. 



4.1 Use of FRE in Chemical Engineering 

The use of fuzzy relational equations (FRE) for the first time 
has been used in the study of flow rates in chemical plants. They 
have only used the concept of linear algebraic equations to 
study this problem and have shown that use of linear equations 
does not always guarantee them with solutions. Thus we are not 
only justified in using fuzzy relational equation but we are 
happy to state by adaptation of FRE we are guaranteed of 
solutions to the problem. We have adapted the fuzzy relational 
equations to the problem of estimation of flow rates in a 
chemical plant, flow rates in a pipe network and use of FRE in a 
3 stage counter current exaction unit [44]. 
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Experimental study of chemical plants is time consuming 
expensive and need intensive labor, researchers and engineers 
prefer only theoretical approach, which is inexpensive and 
effective. Only linear equations have been used to study: (1). A 
typical chemical plant having several inter-linked units (2). 
Flow distribution in a pipe network and (3). A three stage 
counter current extraction unit. Here, we tackle these problems 
in 2 stages. At the first stage we use FRE to obtain a solution. 
This is done by the method of partitioning the matrix as rows. If 
no solution exists by this method we as the second stage adopt 
Fuzzy Neural Networks by giving weightages. We by varying 
the weights arrive at a solution which is very close to the 
predicted value or the difference between the estimated value 
and the predicted value is zero. Thus by using fuzzy approach 
we see that we are guaranteed of a solution which is close to the 
predicted value, unlike the linear algebraic equation in which 
we may get a solution and even granted we get a solution it may 
or may not concur with the predicted value. 

To attain both solution and accuracy we tackle the problems 
using Fuzzy relational equations at the first stage and if no 
solution is possible by this method we adopt neural networks at 
the second stage and arrive at a solution. 

Consider the binary relation P(X, Y), Q(Y, Z) and R(X, Z) 
which are defined on the sets X = {x ; / i e 1} Y = {y ; / j e J} 
and Z{z k / k e K} where we assume that I = N n , J = N r and K = 
N s . Let the membership matrices of P, Q and R be denoted by P 
= [Pij], Q = [qik] and R = [r ik ] respectively, where p tj = P(x ; , yj), 
Qik = Q(yj, Zk) and r ik = R(x ; , z k ) for i e I (= N n ), j e J (= N m ) and 
k g K (= N s ). Entries in P, Q and R are taken from the interval 
[0, 1]. The three matrices constrain each other by the equation 

P o Q = R (1) 

(where o denotes the max-min composition) known as the fuzzy 
relation equation (FRE) which represents the set of equation 

Max pijq jk = r ik (2) 
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for all i e N n , k e N s . If after partitioning the matrix and 
solving the equation (1) yields maximum of q jk < r ik for some 
qjk, then this set of equation has no solution. So at this stage to 
solve the equation 2, we use feed-forward neural networks of 
one layer with n-neurons with m inputs shown in Figure 4.1.1. 

Inputs of the neuron are associated with real numbers Wy 
referred as weights. The linear activation function f is defined 

by 
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Figure: 4.1.1 

The output y ; = f(max WyXj), for all i e N n and j e N m . 
Solution to (1) is obtained by varying the weights Wy so that the 
difference between the predicted value and the calculated value 
is zero. 

FRE to estimate flow rates in a chemical plants 

A typical chemical plant consists of several interlinked units. 
These units act as nodes. The flowsheet is given in Figure 4.1.2. 



81 



D 2 



Di 



M 



Figure: 4.1.2 

An experimental approach would involve measuring the 
nine flow-rates to describe the state of the plant which would 
involve more money and labor. 

While studying this problem in practice researchers have 
has neglected density variations across each stream. The mass 
balance equations across each node at steady state can be 
written as 



p 3 _ p 2 = p 1) 

F 2 - F 4 = F 5 , 

F4 — F7 = F6, 
F2 + Fs = F5, 

Fs = F 9 - F 6 . 



(3) 



Here F ; represents the volumetric flow rate of the i th stream. 
In equation (3) at least four variables have to be specified or 
determined experimentally. 

The remaining five can then be estimated from the equation 
(3), which is generated by applying the principle of 
conservation of mass to each unit. We assume F b F 5 , F 6 and F 9 
are experimentally measured, equation (3) reads with known 
values on the right-hand side as follows: 
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(5) 



where, P, Q and R are explained. Using principle of 
conservation of mass balanced equation we estimate the flow 
rates of the five liquid stream. We in this problem aim to 
minimize the errors between the measured and the predicted 
value. We do this by giving suitable membership grades py e [0, 
1] and estimate the flow rates by using these p^'s in the equation 
3. Now the equation 4 reads as follows: 
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where P = (py), 

Q = (q ik ) = [F 2 F 3 F 4 F 7 F 8 ] ( and 

R = (r ik ) = [F 1 F 5 F 6 F 5 F 9 -F 6 ] t . 
We now apply the partitioning method of solution to 
equation (6). The partitioning of P correspondingly partitions R, 
which is give by a set of give subsets as follows: 
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[0000 p 55 ] 



"F 2 1 *1 

F 4 = F 6 
F 5 F 5 

F F — F 

_8 J r 9 r 6 



Suppose the subsets satisfies the condition max q ik < r ik then it 
has no solution. If it does not satisfy, this condition, then it has a 
final solution. If we have no solution we proceed to the second 
stage of solving the problem using Fuzzy Neural Networks. 

When the FRE has no solution by the partition method, we 
solve these FRE using neural networks. This is done by giving 
weightages of zero elements as and the modified FRE now 
reads as 



F 2 F t 

F 3 F 5 

F 4 = F 6 
F 7 F 5 

^8 ^9 _ ^6 



The linear activation function f defined earlier gives the output 
y ; = f (max Wy Xj) (i e N n ) we calculate max WyXj as follows: 

1. WnXi = 0.02F 2 , Wi 2 x 2 = 0F 2 , W13X3 = 0F 2 W14X4 = 0.045F 2 , 
W 15 x 5 = 0F 2 

yi = f (max J6Nm WyXj) = f (0.02F 2 , 0F 2 , 0.045F 2 , 0F 2 ) 

2. W 2 iXi = 0.04F 3 , W 22 x 2 = 0.045F 3 , W 23 x 3 = 0F 3 , W 24 x 4 = 0.0F 3 , 
Wi 5 x 5 = 0F 3 

y 2 = f (max jeNm W ijXj ) = f (0.04F 3 , 0.045F 3 , 0F 3 , 0.0F 3 , 0F 3 ) 

3. W 3 iX! = O.OF4, W 32 x 2 = 0.085F 4 , W 33 x 3 = 0.15F 4 , W 34 x 4 = 
0.0F 4 W 35 x 5 = 0F 4 

y 3 = f (max J6Nm WyXj) = f (0F 4 , 0.085F 4 , 0.15F 4 , 0F 4 , 0F 4 ) 
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4. W41X1 = 0.0F 7 , W 42 x 2 = 0F 7 , W43X3 = 0.2F 7 , W 44 x 4 = 0.0F 7 , 
W 45 x 5 = 0F 7 

y 4 = f (max jeNm W ijXj ) = f (0F 7 , 0F 7 , 0.2F 7) 0.0F 7) 0F 7 ) 

5. W51X1 = 0.0F 8) W 52 x 2 = 0F 8) W53X3 = 0F 8) W 54 x 4 = 0.45F 8 , 
W 55 x 5 = 0.5F 8 

y 5 = f (max jeNm W ijXj ) = f (0F 8 , 0F 8) 0F 8) 0.45F 8) 0.5F 8 ) 

shown in Figure 4.1.2. Suppose the error does not reach we 
change the weights till the error reaches 0. We continue the 
process again and again until the error reaches to zero. 

Thus to reach the value zero we may have to go on giving 
different weightages (finite number of time) till say s th stage P s o 
Q s whose linear activation function f, makes the predicted value 
to be equal to the calculated value. Thus by this method, we are 
guaranteed of a solution which coincides with the predicted 
value. 



4.2 Fuzzy neural networks to estimate velocity of flow 
distribution in a pipe network 

In flow distribution in a pipe network of a chemical plant, we 
consider liquid entering into a pipe of length T and diameter D 
at a fixed pressure P ; , The flow distributes itself into two pipes 
each of length Ti(T 2 ) and diameter Di(D 2 ) given in Figure 4.2.1. 



Pa, Dl Vi 




Pa, D 2 , V 2 



Figure: 4.2.1 

The linear equation is based on Ohm's law, the drop in 
voltage V across a resistor R is given by the linear relation V = 
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iR (Ohm's law). The hydrodynamic analogue to the mean 
velocity v for laminar flow in a pipe is given by V p = v 
(32liT/D 2 ). This is classical-Poiseulle equation. In flow 
distribution in a pipe network, neglecting pressure losses at the 
junction and assuming the flow is laminar in each pipe, the 
macroscopic momentum balance and the mass balance at the 
junction yields, 

Pi - P a = (32liT/D 2 )v + (32liT 1 D 1 2 )v 1 , 



Pi - P a = (32liT/D 2 )v + 32liT/D 2 2 )v 2 , 
D 2 v = Di 2 vi + v 2 D 2 2 . 



(1) 



Hence P a is the pressure at which the fluid leaves the system at 
the two outlets. The set of three equation in (1) can be solved 
and we estimate v, Vi, v 2 for a fixed (Pi - P a ). The system reads 
as 
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We transform this equation into a fuzzy relation equation. We 
use a similar procedure described earlier and obtain the result by 
fuzzy relation equation. We get max (0.2v, 0.025v, 0.03v), max 
(0.035v, Ovi, 0.04vi), max (0v 2 , 0.04v 2 , 0.045v 2 ) by using neural 
networks for fuzzy relation equation described in [11]. Suppose 
the error does not reach to 0, we change the weights till the error 
reaches 0. We continue the process again and again till the error 
reaches zero. 



4.3 Fuzzy neural networks to estimate three stage counter 
current extraction unit 

Three-stage counter extraction unit is shown in Figure 4.3.1. 
The components A present in phase E (extract) along with a 
nondiffusing substance as being mixture. 
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Stage 1 Stage 2 Stage 3 
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Figure: 4.3.1 

It is extracted into Rby a nondiffusing solvent. The 3 extraction 
stage is given by the three equation. 



E S Y 4 + R S X S = R S X 3 + E S Y 3) 
EsY 3 + RsXi = E s + R S X 2 , 
EsY2 + R s Xq = E s Yi + R s Xi 



(1) 



Yi(Xi) = moles of A, The flow of each stage is denoted by E S (R S ) 
and this constant does not vary between the different stages. The 
assumption of a linear equilibrium relationship for the 
compositions leaving the i th stage equations 



Y ; = KX, 
for i = 1, 2, 3 reads as 



(2) 
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where {Xi, Y h X 2 , Y 2 , X 3 , Y 3 } can be obtained for a given E s , 
R s and K. Since use of linear algebraic equation does not result 
in the closeness of the measured and predicted value, we use 
neural networks for fuzzy relation equations to estimate the 
flow-rates of the stream, moles of the three-stage counter 
extraction unit and velocity of the flow distribution in a pipe 
network. As neural networks is a method to reduce the errors 
between the measured value and the predicted value. This 
allows varying degrees of set membership (weightages) based 
on a membership function defined over the range of value. The 
(weightages) membership function usually varies from to 1. 
We use the similar modified procedure described earlier and get 
result by fuzzy relation equation. We get max (0.2Xi, 0.25Xi, 
0.3Xi, OXi, OXi, OXi), max (0.35Y b 0.4Y b 0Y 1; 0Y 1; 0Y b OYi) 
max (0X 2 , 0X 2 , 0.45X 2 , 0.5X 2 , 0.55X 2 , 0X 2 ), max (0.6Y 2 , 0Y 2 , 
0.65Y 5 , 0.7Y 2 , 0Y 2 , 0Y 2 ) max (0X 3 , 0X 3 , 0X 3 , 0X 3 , 0.75X 3 , 
0.8X 3 ), max (0Y 3 , 0Y 3 , 0.85Y 3 , 0Y 3 , 0.9Y 3 , 0.95Y 3 ) by neural 
networks for fuzzy relation equation. We continue this process 
until the error reaches zero or very near to zero. 

Thus we see that when we replace algebraic linear equations by 
fuzzy methods to the problems described we are not only 
guaranteed of a solution, but our solution is very close to the 
predicted value. 
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Chapter Five 



Mnimzation of \Npste Gas Flow in 
Cheivical Industries 



Chemical Industries and Automobiles are extensively 
contributing to the pollution of environment, Carbon monoxide, 
nitric oxide, ozone, etc., are understood as the some of the 
factors of pollution from chemical industries. The maintenance 
of clean and healthy atmosphere makes it necessary to keep the 
pollution under control which is caused by combustion waste 
gas. The authors have suggested theory to control waste gas 
pollution in environment by oil refinery using fuzzy linear 
programming. To the best of our knowledge the authors [43]are 
the first one to apply fuzzy linear programming to control or 
minimize waste gas in oil refinery. 

An oil refinery consists of several inter linked units. These 
units act as production units, refinery units and compressors 
parts. These refinery units consume high-purity gas production 
units. But the gas production units produce high-purity gas 
along with a low purity gas. This low purity gas goes as a waste 
gas flow and this waste gas released in the atmosphere causes 
pollution in the environment. But in the oil refinery the quantity 
of this waste gas flow is an uncertainty varying with time and 
quality of chemicals used in the oil refinery. Since a complete 
eradication of waste gas in atmosphere cannot be made; here 
one aims to minimize the waste gas flow so that pollution in 
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environment can be reduced to some extent. Generally waste 
gas flow is determined by linear programming method. In the 
study of minimizing the waste gas flow, some times the current 
state of the refinery may already be sufficiently close to the 
optimum. To over come this situation we adopt fuzzy linear 
programming method. 
The fuzzy linear programming is defined by 

Maximize z = ex 

Such that Ax < b 
x<0 
where the coefficients A, b and c are fuzzy numbers, the 
constraints may be considered as fuzzy inequalities with 
variables x and z. We use fuzzy linear programming to 
determine uncertainty of waste gas flow in oil refinery which 
pollutes the environment. 

Oil that comes from the ground is called "Crude oil". 
Before one can use it, oil has to be purified at a factory called a 
"refinery", so as to convert into a fuel or a product for use. The 
refineries are high-tech factories, they turn crude oil into useful 
energy products. During the process of purification of crude oil 
in an oil refinery a large amount of waste gas is emitted to 
atmosphere which is dangerous to human life, wildlife and plant 
life. The pollutants can affect the health in various ways, by 
causing diseases such as bronchitis or asthma, contributing to 
cancer or birth defects or perhaps by damaging the body's 
immune system which makes people more susceptible to a 
variety of other health risks. Mainly, this waste gas affects 
Ozone Layer. Ozone (or Ozone Layer) is 10-50 km above the 
surface of earth. Ozone provides a critical barrier to solar 
ultraviolet radiation, and protection from skin cancers, cataracts, 
and serious ecological disruption. Further sulfur dioxide and 
nitrogen oxide combine with water in the air to form sulfuric 
acid and nitric acid respectively, causing acid rain. It has been 
estimated that emission of 70 percentage of sulfur dioxide and 
nitrogen oxide are from chemical industries. 

We cannot stop this process of oil refinery, since oil and 
natural gas are the main sources of energy. We cannot close 
down all oil refineries, but we only can try to control the amount 
of pollution to a possible degree. In this paper, the authors use 
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fuzzy linear programming to reduce the waste gas from oil 
refinery. The authors describe the knowledge based system 
(KBS) that is designed and incorporate it in this paper to 
generate an on-line advice for operators regarding the proper 
distribution of gas resources in an oil refinery. In this system, 
there are many different sources of uncertainty including 
modeling errors, operating cost, and different opinions of 
experts on operating strategy. The KBS consists of sub- 
functions, like first sub-functions, second sub-functions, etc. 
Each and every sub-functions are discussed relative to certain 
specific problems. 

For example: The first sub-function is mainly adopted to the 
compressor parts in the oil refineries. Till date they were using 
stochastic programming, flexibility analysis and process design 
problems for linear or non-linear problem to compressor parts in 
oil refinery. Here we adopt the sub function to study the proper 
distribution of gas resources in an oil refinery and also use fuzzy 
linear programming (FLP) to minimize the waste gas flow. By 
the term proper distribution of gas we include the study of both 
the production of high-purity gas as well as the amount of waste 
gas flow which causes pollution in environment. 

In 1965, Lofti Zadeh [115, 116] wrote his famous paper 
formally defining multi-valued, or "fuzzy" set theory. He 
extended traditional set theory by changing the two-values 
indicator functions i.e., 0, 1 or the crisp function into a multi- 
valued membership function. The membership function assigns 
a "grade of membership" ranging from to 1 to each object in 
the fuzzy set. Zadeh formally defined fuzzy sets, their 
properties, and various properties of algebraic fuzzy sets. He 
introduced the concept of linguistic variables which have values 
that are linguistic in nature (i.e. pollution by waste gas = {small 
pollution, high pollution, very high pollution}). 

Fuzzy Linear Programming (FLP): FLP problems with 
fuzzy coefficients and fuzzy inequality relations as a multiple 
fuzzy reasoning scheme, where the past happening of the 
scheme correspond to the constraints of thee FLP problem. We 
assign facts (real data from industries) of the scheme, as the 
objective of the FLP problem. Then the solution process 
consists of two steps. In the fist step, for every decision 
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variable, we compute the (fuzzy) value of the objective function 
via constraints and facts/objectives. At the second step an 
optimal solution to FLP problem is obtained at any point, which 
produces a maximal element to the set of objective functions (in 
the sense of the given inequality relation). 

The Fuzzy Linear Programming (FLP) problem application 
is designed to offer advice to operating personnel regarding the 
distribution of Gas within an oil refinery (Described in Figure 
5.1) in a way which would minimize the waste gas in 
environment there by reduce the atmospheric pollution . 

GPUI, GPU2 and GPU3 are the gas production units and 
GGG consumes high purity gas and vents low purity gas. Gas 
from these production units are sent to some oil refinery units, 
like sulfur, methanol, etc. Any additional gas needs in the oil 
refinery must be met by the gas production unit GPU3. 

The pressure swing adsorption unit (PSA) separates the 
GPU2 gas into a high purity product stream and a low purity tail 
stream (described in the Figure 5.1). Ci, C 2 , C 3 , C 4 , C 5 , are 
compressors. The flow lines that dead -end is an arrow 
represent vent to flare or fuel gas. This is the wasted gas that is 
to be minimized. Also we want to minimize the letdown flow 
from the high purity to the low purity header 
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Figure: 5.1 
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FLP is a method of accounting for uncertainty is used by the 
authors for proper distribution of gas resources, so as to 
minimize the waste gas flow in atmosphere. FLP allows varying 
degrees of set membership based on a membership function 
defined over a range of values. The membership function 
usually varies from to 1. FLP allow the representation of 
many different sources of uncertainty in the oil refinery. These 
sources may (or) may not be probabilistic in nature. The 
uncertainty is represented by membership functions describing 
the parameters in the optimization model. A solution is found 
that either maximizes a given feasibility measure and 
maximizes the wastage of gas flow. FLP is used here to 
characterize the neighborhood of solutions that defines the 
boundaries of acceptable operating states. 

Fuzzy Linear Programming (FLP) can be stated as; 



max imize z = ex 

s.t Ax < b 

x>0 



•••(*) 



The coefficients A, b and c are fuzzy numbers, the 
constraints may be considered as fuzzy inequalities. The 
decision space is defined by the constraints with c, x e N, b e 
R m and A e R m , where N, R m , and R mxn are reals. 

The optimization model chosen by the knowledge based 
system (KBS) is determined online and is dependent on the 
refinery units. This optimization method is to reduce the amount 
of waste gas in pollution. 
We aim to 

1. The gas (GCG 2 ) vent should be minimized. 

2. The let down flow should be minimized and 

3. The make up gas produced by the as production unit 
(GPU3) should be minimized. 
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Generally the waste gas emitted by the above three ways 
pollute the environment. The objective function can be 
expressed as the sum of the individual gas waste flows. The 
constrains are given by some physical limitations as well as 
operator entries that describe minimum and maximum desired 
flows. 

The obtained or calculated resultant values of the decision 
variables are interpreted as changes in the pressure swing 
adsorption feed, and the rate that gas is imported to CGG and 
gas production unit (GPU3). But in the optimization model 
there is uncertainty associated with amount of waste gas from 
oil refinery, and also some times the current state of the refinery 
may already be sufficiently close to the optimum. 

For example to illustrate the problem, if the fuzzy 
constraints x u the objects are taken along the x-axis are shown 
in the figures 5.2 and 5.3, which represent the expression. 




Figure: 5.2 




Figure: 5.3 



Xi < 8 (with tolerance p = 2) 



(1) 



The membership function p are taken along the y-axis i.e. 
p(xx) lies in [0, 1] this can be interpreted as the confidence with 
which this constraint is satisfied (0 for low and 1 for high). The 
fuzzy inequality constraints can be redefined in terms of their a- 
cuts. 

{S a / a 8 [0, 1]}, where S a = {y / (p (y) > a)}. 

The parameter a is used to turn fuzzy inequalities into crisp 
inequalities. So we can rewrite equation (1) 

xi < 6 + 2 (2) (1- a) 
X! < 6 + 4 (1 - a) 
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where a 8 [0, 1] expressed in terms of a in this way the fuzzy 
linear programming problem can be solved parametrically. The 
solution is a function on a 

x* = f(a) (2) 

with the optimal value of the objective function determined by 
substitution in equation (1). 

z* = ex* = g(a). (3) 

This is used to characterize the objective function. The result 
covers all possible solutions to the optimization problem for any 
point in the uncertain interval of the constraints. 

The a-cuts of the fuzzy set describes the region of feasible 
solutions in figures 5.2 and 5.3. The extremes (a = and a = 1) 
are associated with the minimum and maximum values of x 
respectively. The given equation (2) can also be found this, is 
used to characterize the objective function. The result covers all 
possible solutions to the optimization problem for any point in 
the uncertain interval of the constraints. 

Fuzzy Membership Function to Describe Uncertainty: The 
feasibility of any decision (p D ) is given by the intersection of 
the fuzzy set describing the objective and the constraints. 

p D (x) = p z (x) A p N (x) 

where A represents the minimum operator, that is the usual 
operation for fuzzy set intersection. The value of p N can be 
easily found by intersecting the membership values for each of 
the constraints. 

p N (x) = p 1 (x)Ap 2 ( X )A...A Um ( x ). 

The membership functions for the objective (p z ) however is 
not obvious z is defined in (2). Often, predetermined aspiration 
target values are used to define this function. Since reasonable 
values of this kind may not be available, the solution to the FLP 
equation (3) is used to characterize this function. 
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Pz(x) = 



1 

z(x)-b(l) 

b(0)-b(l) 





if z(x)>b(0) 

if b(l) < z(x) < b(0) 

if z(x) < b(0). 



(5) 



The result is that the confidence value increases as the value of 
the objective value increases. This is reasonable because the 
goal is to maximize this function the limits on the function 
defined by reasonable value is obtained by extremes of the 
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Figure: 5.4 



object 



objective value. 



These are the results generated by the fuzzy linear 
programming. Since both u N and u z have been characterized, 
now our goal is to describe the appropriateness of any operation 
state. Given any operating x, the feasibility can be specified 
based on the objective value, the constraints and the estimated 
uncertainty is got using equation (4). The value of p D are shown 
as the intersection of the two membership functions. 

Defining the decision region based on the intersection we 
describe the variables and constraints of our problem. The 
variable Xi represents the amount of gas fed to pressure swing 
adsorption from the gas production unit. The variable x 2 
represents the amount of gas production that is sent to CGG. 
This problem can be represented according to equation (*). The 
constraints on the problem are subjected to some degree of 
uncertainty often some violation of the constraints within this 
range of uncertainty is tolerable. This problem can be 
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represented according to equation (*). Using the given refinery 
data from the chemical plant. 
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Figure: 5.5 



object 



c = [-0.544 3] 



A = 



1 


0~ 





1 


-0.544 


1 



b = 



33.652 
23.050 
4.743 



Using equation (*) we get 

Zc = - 0.544 Xi + 3x 2 it represents gas waste flow. The gas 

waste flow is represented by the following three equations: 



l. 

ii. 



m. 



Xi + 0x 2 < 33.652 is the total dead - end waste flow gas. 

Oxi + x 2 = 23.050 is the total (GCG2) gas consuming 

gas - treaters waste flow gas. 

- 0544 X! + x 2 < 4.743 is the total let-down waste flow 

gas. 



All flow rates are in million standard cubic feet per day. (i.e. 1 
MMSCFD = 0.3277 m 3 /s at STP). The value used for may be 
considered to be desired from operator experts opinion. The 
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third constraint represents the minimum let-down flow receiving 
to keep valve from sticking. The value to this limit cannot be 
given an exact value, therefore a certain degree of violation may 
be tolerable. The other constraints may be subject to some 
uncertainty as well as they represent the maximum allowable 
values for Xi and x 2 . In this problem we are going to express all 
constrains in terms of a, a, e [0, 1]. We have to chose a value of 
tolerance on the third constraint as p 3 = 0.1, then this constraint 
is represented parametrically as 

a 3 x < (b 3 - p 3 ) + 2p 3 (1 - a). 

For example, if we use crisp optimization problem with the 
tolerance value p = 0.1 we obtain the following result: 




22.5 



32 



33 33.467 33.5 X: 

Figure: 5.6 



where Xi represents the amount of gas fed to PSA from gas 
production unit which is taken along the x axis, and x 2 amount 
of gas sent to CGG which is taken along the y axis, 
we get xi = 33.469, when x 2 = 23.050 

[33.469" 

[23.050_ 

z = 50.941. Finally we compare this result with our fuzzy linear 
programming method. 

We replace two valued indicator function method by fuzzy 
linear programming. 
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Fuzzy Linear Programming is used now to maximize the 
objective function as well as minimize the uncertainty (waste 
flow gas). For that all of the constraints are expressed in terms 
of a, a, e [0, 1]. 

a 3 x < (b 3 - p 3 ) + 2p 3 (1 - a), a e [0, 1] 
where a 3 is the third row in the matrix A. i.e. = 0.544xi + x 2 < 
4.843 - 0.2 a, when the tolerance p 3 = 0.3, we fix the value of a 
8 [0.9,1], when the tolerance p 3 = 0.1, we see a 8 [0.300, 0.600]. 




22.5 



33.469 



Figure: 5.7 

where xi represents the amount of gas fed to PSA from gas 
production unit which is taken along the x axis, and x 2 amount 
of gas that is sent to CGG which is taken along the y axis, 

When x 2 = 23.050 and a = 0.0, we get Xi = 33.469. 

When x 2 = 23.050 and a = 0.4, we get x l = 33.616 
The set (u z ) is defined in equation 5. Fuzzy Linear 
Programming solution is 

[33.469" 
x* = f(a) = 

[23.050 

this value is recommended as there is no changes in the 
operating policy. 

So we have to chose the value for a as 0.6 for the tolerance 
p 3 = 0.1, we get the following graph where Xi represents the 
amount of gas fed to PSA from gas production unit which is 
taken along the x axis, and x 2 amount of gas that sent to CGG 
which is taken along the y axis, 
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22.5 



Figure: 5.8 

when x 2 = 23.050 and a = 0.0 we get x a = 33.469 

when x 2 = 23.050 and a = 0.6 we get Xj = 33.689. 

The operating region 

"33.689" 
x* = f (0.6) = 

[23.050 

Now if the tolerance on the third constraint is increased to p 3 = 
0.2. This results is the region shown in the following graph. As 
expected the region has increased to allow a larger range of 
operating states. 

when x 2 = 23.050 and a = 0.0 we get x l = 33.285 
when x 2 = 23.050 and a = 0.9 we get x l = 33.947. 
The operating region is 




31.5^ 



33.285 33.5 



Figure: 5.9 
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where Xi represents the amount of gas fed to PSA from gas 
production unit which is taken along the x axis, and x 2 amount 
of gas that is sent to CGG which is taken along the y axis. 



x* = f (0.9) = 



33.947 
23.050 



The fuzzy linear programming solution is 

"33.285 
t* = f (a) 

23.050 



z* = 51.043. 

Finally we have to take a 8 [0.9, 1.00]. 

Choose a = and when the tolerance p 3 = 0.3 we get the 

following graph when x 2 = 23.050 we get Xi = 33.101. 




33.101 33.5 Xi 



Figure: 5.10 



where x 1 represents the amount of gas fed to PSA from gas 
production unit; and x 2 amount of gas that is sent to CGG. 

When a = 1 and x 2 = 23.050 we get x x = 34.204. The 
operating region is 



f (1-0) 



33.204 
23.050 
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The fuzzy linear programming solutions are 
x* = f (a) 



33.101 
23.050 



The fuzzy linear programming solutions are 
z* = g (a) = 51.143. 

We chose maximum value from the Fuzzy Linear Programming 
method i.e. z* = 51.143. 

Thus when we work by giving varying membership 
functions and use fuzzy linear programming we see that we get 
the minimized waste gas flow value as 33.101 in contrast to 
33.464 measured in million standard cubic feet per day and the 
maximum gas waste flow of system is determined to be 51.143 
in contrast to their result of 50.941 measured in million standard 
cubic feet per day. Since the difference we have obtained is 
certainly significant, this study when applied to any oil refinery 
will minimize the waste gas flow to atmosphere considerably 
and reduce the pollution. 
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Chapter Six 



Use of Neutrosophic Relational 
Equations in Chemical Engineering 



This chapter has 2 sections. Section one gives introduction to 
Neutrosophic Relational Equations (NRE) and section two gives 
use of NRE in chemical Engineering. 



6.1 Introduction to Neutrosophic Relation and their 
properties 

In this section we introduce the notion of neutrosophic relational 
equations and fuzzy neutrosophic relational equations and 
analyze and apply them to real-world problems, which are 
abundant with the concept of indeterminacy. We also mention 
that most of the unsupervised data also involve at least to certain 
degrees the notion of indeterminacy. 

Throughout this section by a neutrosophic matrix we mean 
a matrix whose entries are from the set N = [0, 1] u I and by a 
fuzzy neutrosophic matrix we mean a matrix whose entries are 
from N* = [0, l]u{nl/ne (0,1]}. 

Now we proceed on to define binary neutrosophic relations 
and binary neutrosophic fuzzy relation. 
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A binary neutrosophic relation R N (x, y) may assign to each 
element of X two or more elements of Y or the indeterminate I. 
Some basic operations on functions such as the inverse and 
composition are applicable to binary relations as well. Given a 
neutrosophic relation R N (X, Y) its domain is a neutrosophic set 
onXul domain R whose membership function is defined by 

domR(x) = max R N (x, y) 

yeXuI 

for each xeXuI. 

That is each element of set X u J belongs to the domain of 
R to the degree equal to the strength of its strongest relation to 
any member of set YuJ. The degree may be an indeterminate I 
also. Thus this is one of the marked difference between the 
binary fuzzy relation and the binary neutrosophic relation. The 
range of R N (X,Y) is a neutrosophic relation on Y, ran R whose 
membership is defined by 

ran R(y) = max R N (x, y) 

xeX 

for each ye Y, that is the strength of the strongest relation that 
each element of Y has to an element of X is equal to the degree 
of that element's membership in the range of R or it can be an 
indeterminate I. 

The height of a neutrosophic relation R N (x, y) is a number 
h(R) or an indeterminate I defined by 

h N (R) = max max R N (x, y). 

yeYuI xeXuI 

That is h N (R) is the largest membership grade attained by any 
pair (x, y) in R or the indeterminate I. 

A convenient representation of the neutrosophic binary 
relation Rn(X, Y) are membership matrices R = [y xy ] where y xy 
e R N (x, y). 

Another useful representation of a binary neutrosophic 
relation is a neutrosophic sagittal diagram. Each of the sets X, Y 
represented by a set of nodes in the diagram, nodes 
corresponding to one set are clearly distinguished from nodes 
representing the other set. Elements of X' x Y' with non-zero 
membership grades in Rn(X, Y) are represented in the diagram 
by lines connecting the respective nodes. These lines are labeled 
with the values of the membership grades. 
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An example of the neutrosophic sagittal diagram is a binary 
neutrosophic relation Rn(X, Y) together with the membership 
neutrosophic matrix is given below. 



yi y2 ys y4 

10 0.5 

0.3 0.4 

1 0.2 

10 

0.5 0.7 



The inverse of a neutrosophic relation Rn(X, Y) = R(x, y) 
for all x e X and all y e Y. A neutrosophic membership matrix 
R _1 = [r^ 1 ] representing R^Y, X) is the transpose of the 

matrix R for Rn(X, Y) which means that the rows of R" 1 equal 
the columns of R and the columns of R" 1 equal rows of R. 
Clearly (R 1 )" 1 = R for any binary neutrosophic relation. 

Consider any two binary neutrosophic relations P N (X, Y) 
and Q N (Y, Z) with a common set Y. The standard composition 
of these relations which is denoted by P N (X, Y) • Qn(Y, Z) 
produces a binary neutrosophic relation Rn(X, Z) on X x Z 
defined by 
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R N (x, z) = [P • Q] N (x, z) = max min[P N (x, y), Q N (x, y)] 

ysY 



for all xe X and all z e Z. 

This composition which is based on the standard t N -norm 
and t N -co-norm, is often referred to as the max-min 
composition. It can be easily verified that even in the case of 
binary neutrosophic relations 

[P N (X, Y) • Q N (Y, Z)] 1 

= Q^(Z, Y) • P^(Y, X). [P N (X, Y) • Q N (Y, Z)] • R N (Z, W) 

= P N (X, Y) • [Q N (Y, Z) • R N (Z, W)], 

that is, the standard (or max-min) composition is associative 
and its inverse is equal to the reverse composition of the inverse 
relation. However, the standard composition is not 
commutative, because Q N (Y, Z) • P N (X, Y) is not well defined 
when X * Z. Even if X = Z and Q N (Y, Z) ° P N (X, Y) are well 
defined still we can have P N (X, Y) ° Q (Y, Z) * Q (Y, Z) ° P 
(X, Y). 

Compositions of binary neutrosophic relation can the 
performed conveniently in terms of membership matrices of the 
relations. Let P = [p ik ], Q = [q kj ] and R = [ry] be membership 
matrices of binary relations such that R = P ° Q. We write this 
using matrix notation 

fry] = [pik] o [q kj ] 
where ry = maxmin (pi k , qkj). 

k 

A similar operation on two binary relations, which differs 
from the composition in that it yields triples instead of pairs, is 
known as the relational join. For neutrosophic relation P N (X, Y) 
and Q N (Y, Z) the relational join P * Q corresponding to the 
neutrosophic standard max-min composition is a ternary relation 
R N (X, Y, Z) defined by R N (x, y, z) = [P * Q] N (x, y, z) = min 
[P N (x, y), Q N (y, z)] for each x e X, y e Y and z e Z. 

This is illustrated by the following Figure 6.1.2. 
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In addition to defining a neutrosophic binary relation there 
exists between two different sets, it is also possible to define 
neutrosophic binary relation among the elements of a single set 
X. A neutrosophic binary relation of this type is denoted by 
R N (X, X) or R N (X 2 ) and is a subset of X x X = X 2 . 

These relations are often referred to as neutrosophic 
directed graphs or neutrosophic digraphs. [42] 

Neutrosophic binary relations R N (X, X) can be expressed 
by the same forms as general neutrosophic binary relations. 
However they can be conveniently expressed in terms of simple 
diagrams with the following properties. 

I. Each element of the set X is represented by a single 
node in the diagram. 

II. Directed connections between nodes indicate pairs of 
elements of X for which the grade of membership in 
R is non zero or indeterminate. 

III. Each connection in the diagram is labeled by the 
actual membership grade of the corresponding pair in 
R or in indeterminacy of the relationship between 
those pairs. 

The neutrosophic membership matrix and the neutrosophic 
sagittal diagram is as follows for any set X = {a, b, c, d, e}. 
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Neutrosophic membership matrix for x is given above and the 
neutrosophic sagittal diagram is given below. 



a 
b 
c 
d 


8sS 


a 

b 
c 

d 


e 


■^•■* < "^. ^v^ 


e 


• .2 ^» 

Figure: 4.1.3 



Neutrosophic diagram or graph is left for the reader as an 
exercise. 

The notion of reflexivity, symmetry and transitivity can be 
extended for neutrosophic relations R N (X, Y) by defining them 
in terms of the membership functions or indeterminacy relation. 

Thus R N (X, X) is reflexive if and only ifR N (x, x) = 1 for all x £ 
X. If this is not the case for some x e X the relation is 
irreflexive. 

A weaker form of reflexivity, if for no x in X; Rn(x, x) = 1 then 
we call the relation to be anti-reflexive referred to as e- 
reflexivity, is sometimes defined by requiring that 
R N (x, x) > e where < e <1. 

A fuzzy relation is symmetric if and only if 

R N (x, y) = R N (y, x) for all x, y, e X. 
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Whenever this relation is not true for some x, y e X the relation 
is called asymmetric. Furthermore when R N (x, y) > and R N (y, 
x) > implies that x = y for all x, y e X the relation R is called 
anti-symmetric. 

A fuzzy relation R N (X, X) is transitive (or more specifically 
max-min transitive) if 

R N (x, z) > max min [R N (x, y), R N (y, z)] 

is satisfied for each pair (x, z) e X 2 . A relation failing to satisfy 
the above inequality for some members of X is called non- 
transitive and if R N (x, x) < max min [Rn(x, y), Ri^y, z)] for all 

ysY 

(x, x) e X 2 , then the relation is called anti-transitive 

Given a relation Rn(X, X) its transitive closure R NT (x, X) 
can be analyzed in the following way. 

The transitive closure on a crisp relation R N (X, X) is 
defined as the relation that is transitive, contains 

R N (X, X) < max min [R N (x, y) R N (y, z)] 

yeF 

for all (x, x) e X 2 , then the relation is called anti- transitive. 
Given a relation R N (x, x) its transitive closure R NT (X, X) can 
be analyzed in the following way. 

The transitive closure on a crisp relation R N (X, X) is 
defined as the relation that is transitive, contains R N and has the 
fewest possible members. For neutrosophic relations the last 
requirement is generalized such that the elements of transitive 
closure have the smallest possible membership grades, that still 
allow the first two requirements to be met. 

Given a relation R N (X, X) its transitive closure R NT (X, X) 

can be determined by a simple algorithm. 

Now we proceed on to define the notion of neutrosophic 
equivalence relation. 

Definition 4.1.1: A crisp neutrosophic relation R N (X, X) that 
is reflexive, symmetric and transitive is called a neutrosophic 
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equivalence relation. For each element x in X, we can define a 
crisp neutrosophic set A x which contains all the elements of X 
that are related to x by the neutrosophic equivalence relation. 

Formally A x = [ y | (x, y) € R N (X, X)}. A x is clearly a subset of 
X. The element x is itself contained in A x , due to the reflexivity 
ofR because R is transitive and symmetric each member of A x is 
related to all other members of A x . Further no member of A x is 
related to any element of X not included in A x . This set A x is 
clearly referred to as a neutrosophic equivalence class ofR N (X, 
x) with respect to x. The members of each neutrosophic 
equivalence class can be considered neutrosophic equivalent to 
each other and only to each other under the relation R. 

But here it is pertinent to mention that in some X; (a, b) may 
not be related at all to be more precise there may be an element 
a e X which is such that its relation with several or some 
elements in X \ {a} is an indeterminate. The elements which 
cannot determine its relation with other elements will be put in 
as separate set. 

A neutrosophic binary relation that is reflexive, symmetric 
and transitive is known as a neutrosophic equivalence relation. 

Now we proceed on to define Neutrosophic intersections 
neutrosophic t - norms (t N - norms) 

Let A and B be any two neutrosophic sets, the intersection 
of A and B is specified in general by a neutrosophic binary 
operation on the set N = [0, 1] u I, that is a function of the form 

i N :NxN^N. 
For each element x of the universal set, this function takes as its 
argument the pair consisting of the elements membership grades 
in set A and in set B, and yield the membership grade of the 
element in the set constituting the intersection of A and B. Thus, 

(A n B) (x) = i N [A(x), B(x)] for all x e X. 

In order for the function i N of this form to qualify as a fuzzy 
intersection, it must possess appropriate properties, which 
ensure that neutrosophic sets produced by i N are intuitively 



110 



acceptable as meaningful fuzzy intersections of any given pair 
of neutrosophic sets. It turns out that functions known as t N - 
norms, will be introduced and analyzed in this section. In fact 
the class of t N - norms is now accepted as equivalent to the class 
of neutrosophic fuzzy intersection. We will use the terms t N - 
norms and neutrosophic intersections inter changeably. 

Given a t N - norm, i N and neutrosophic sets A and B we 
have to apply: 

(A nB) (x) = i N [A (x) , B (x)] 

for each x e X, to determine the intersection of A and B based 
upon i N . 

However the function i N is totally independent of x, it 
depends only on the values A (x) and B(x). Thus we may ignore 
x and assume that the arguments of i N are arbitrary numbers 
a, b e [0, 1] u I = N in the following examination of formal 
properties of t N -norm. 

A neutrosophic intersection/ t N -norm i N is a binary operation 
on the unit interval that satisfies at least the following axioms 
for all a, b, c, d e N = [0, 1] u I. 



In 


i N (a, 1) = a 


2 n 


i N (a, I) = I 


3 N 


b < d implies 




i N (a, b) < i N (a, d) 


4 N 


i N (a, b) = i N (b, a) 


5n 


i N (a, i N (b, d)) = i N (a, b), d). 



We call the conditions 1 N to 5 N as the axiomatic skeleton for 
neutrosophic intersections / t N - norms. Clearly i N is a 
continuous function on N \ I and i N (a, a) < a Va e N \ I 

In (II) = L 

If ai < a 2 and bi < b 2 implies i N (a b bi) < i N (a 2 , b 2 ). Several 
properties in this direction can be derived as in case of t-norms. 

The following are some examples of t N -norms 
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1. i N (a, b) = min (a, b) 

i N (a, f) = min (a, T) = I will be called as standard 
neutrosophic intersection. 

2. i N (a, b) = ab for a, b e N \ I and i N (a, b) = I for a, b e 
N where a = I or b = I will be called as the neutrosophic 
algebraic product. 

3. Bounded neutrosophic difference. 

i N (a, b) = max (0, a + b - 1) for a, b e I 
i N (a, 7) = I is yet another example of t N - norm. 
1. Drastic neutrosophic intersection 



2. 



i N (a, b) 



a when b = 1 

b when a = 1 

I when a = I 
orb = I 
or a = b = I 

otherwise 



As I is an indeterminate adjoined in t N - norms. It is not easy to 
give then the graphs of neutrosophic intersections. Here also we 
leave the analysis and study of these t N - norms (i.e. 
neutrosophic intersections) to the reader. 

The notion of neutrosophic unions closely parallels that of 
neutrosophic intersections. Like neutrosophic intersection the 
general neutrosophic union of two neutrosophic sets A and B is 
specified by a function 

)i N :NxN->N where N = [0, 1] u 7. 

The argument of this function is the pair consisting of the 
membership grade of some element x in the neutrosophic set A 
and the membership grade of that some element in the 
neutrosophic set B, (here by membership grade we mean not 
only the membership grade in the unit interval [0, 1] but also the 
indeterminacy of the membership). The function returns the 
membership grade of the element in the set A u B. 

Thus (A u B) (x) = u N [A (x), B(x)] for all x e X. 
Properties that a function u N must satisfy to be initiatively 
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acceptable as neutrosophic union are exactly the same as 
properties of functions that are known. Thus neutrosophic union 
will be called as neutrosophic t-co-norm; denoted by t N - co- 
norm. 

A neutrosophic union / t N - co-norm u N is a binary 
operation on the unit interval that satisfies at least the following 
conditions for all a, b, c, d e N = [0, 1] u I 



Cx 


u N (a, 0) = a 


c 2 


u N (a, 7) = I 


c 3 


b < d implies 




u N (a, b) < u N (a, d) 


c 4 


p N (a, b) = u N (b, a) 


c 5 


p N (a, u N (b, d)) 




= u N (u N (a, b), d). 



Since the above set of conditions are essentially neutrosophic 
unions we call it the axiomatic skeleton for neutrosophic unions 
/ t N -co-norms. The addition requirements for neutrosophic 
unions are 

i. u N is a continuous functions on N \ {/} 

ii. u N (a, a) > a. 

iii. ai < a 2 and bi < b 2 implies u N (ai. bi) < u N (a 2 , b 2 ); 

a 1; a 2 , bi, b 2 e N \ {/} 

We give some basic neutrosophic unions. 
Let u N : [0, 1] x [0, 1] -> [0, 1] 

p N (a, b) = max (a, b) 

p N (a, 7) = I is called as the standard 

neutrosophic union. 
p N (a, b) = a + b - ab and 
u N (a, 1) = I . 

This function will be called as the neutrosophic algebraic sum. 

p N (a, b) = min (1, a + b) and u N (a, I) = I 
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will be called as the neutrosophic bounded sum. We define the 
notion of neutrosophic drastic unions 



Un (a, b) 



a when b = 
b when a = 
I when a = I 
orb = I 
1 otherwise. 



Now we proceed on to define the notion of neutrosophic 
Aggregation operators. Neutrosophic aggregation operators on 
neutrosophic sets are operations by which several neutrosophic 
sets are combined in a desirable way to produce a single 
neutrosophic set. 

Any neutrosophic aggregation operation on n neutrosophic 
sets (n > 2) is defined by a function h N : N n — > N where N = [0, 
l]ul and N n = Nx...xN when applied to neutrosophic sets 

1 v ' 

n- times 

Ai, A 2 ,..., A n defined on X the function h N produces an 
aggregate neutrosophic set A by operating on the membership 
grades of these sets for each x e X (Here also by the term 
membership grades we mean not only the membership grades 
from the unit interval [0, 1] but also the indeterminacy I for 
some x e X are included). Thus 

A N (x) = h N (A! (x), A 2 (x),. . ., A„(x)) 
for each x e X. 



6.2 Use of NRE in chemical engineering 

The use of FRE for the first line has been used in the study of 
flow rates in chemical plants. In this study we are only 
guaranteed of a solution but when we use NRE in study of flow 
rates in the chemical plants we are also made to understand that 
certain flow rates are indeterminate s depending on the leakage, 
chemical reactions and the new effect due to chemical reactions 
which may change due to change in the density/ viscosity of the 
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fluid under study their by changing the flow rates while 
analyzing as a mathematical model. So in the study of flow rates 
in chemical plants some indeterminacy are also related with it. 
FRE has its own limitation for it cannot involve in its analysis 
the indeterminacy factor. 

We have given analysis in chapter 2 using FRE. Now we 
suggest the use of NRE and bring out its importance in the 
determination of flow rates in chemical plants. 

Consider the binary neutrosophic relations P N (X, Y) Q N (Y, 
Z) and R (X, Z) which are defined on the sets X, Y and Z. Let 
the membership matrices of P, Q and R be denoted by P = [py], 
Q = [q jk ] and R = [ry] respectively where py = P(x ; , yj), q jk = Q 
(yj, r k ) and r ik = R (x i( z k ) for ie J = N n , j e J = N m and k e K = N k 
entries of P, Q and R are taken for the interval [0 1] x FN. The 
three neutrosophic matrices constrain each other by the equation 



PoQ = R 



(1) 



where 'o' denotes the max-min composition (1) known as the 
Neutrosophic Relational Equation (NRE) which represents the 
set of equation 

max py q jk = r ik . (2) 



For all i e N n and k e N s . After partitioning the matrix and 
solving the equation (1) yields maximum of q jk < r ik for some 
q jk , then this set of equation has no solution so to solve equation 
(2) we invent and redefine a feed - forward neural networks of 
one layer with n-neurons with m inputs. The inputs are 
associated with Wy called weights, which may be real, or 
indeterminates from RI. The neutrosophic activation function f N 
is defined by 

fO if a < 



f N (a) = 



a if a e [01] 
lif a>l 
al if a e FN 
I if al > I 
if in al, a < 0. 



115 



The out put y ; = f N (max wy Xj). Now the NRE is used to 
estimate the flow rates in a chemical plant. In places where the 
indeterminacy is involved the expert can be very careful and use 
methods to overcome indeterminacy by adopting more and 
more constraints which have not been given proper 
representation and their by finding means to eliminate the 
indeterminacy involved in the weights. 

In case of impossibility to eliminate these indeterminacy 
one can use the maximum caution in dealing with these values 
which are indeterminate s so that all types of economic and time 
loss can be met with great care. In the flow rate problem the use 
of NRE mainly predicts the presence of the indeterminacy 
which can be minimized using f N ; where by all other in- 
descripancies are given due representation. 

We suggest the use of NRE for when flow rates are 
concerned in any chemical plant the due weightage must be 
given the quality of chemicals or raw materials which in many 
cases are not up to expectations, leakage of pipe, the viscosity or 
density after chemical reaction time factor, which is related with 
time temperature and pressure for which usually due 
representations, is not given only ideal conditions are assumed. 
Thus use of NRE may prevent accident, economic loss and 
other conditions and so on. 
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